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Abstract. 

Using the Gaussian+Sch echter composite luminosity functions measured from the ESO-Sculptor Survey 
Ide Lapparent et al.ll2003T) . and assuming that these functions do not evolve with redshift out to z ~ 1, we 
obtain evidence for evolution in the late spectral class containing late-type Spiral (Sc+Sd) and dwarf Irregular 
(di) galaxies. Fhere are indications that the Sc+Sd galaxies are the evolving population, but we cannot exclude 
that the di galaxies also undergo some evolution. Fhis evolution is detected as an increase of the Sc+Sd+Im 
galaxy density which can be modeled as either n(z) tx 1 + 3(z — 0.15) or n{z) oc (1 + z) 2 using the currently 
favored cosmological parameters Q m = 0.3 and SIa = 0.7; the uncertainty in the linear and power-law evolution 
rates is of order of unity. For fl m = 1.0 and Qa = 0.0, the linear and power- law evolution rates are ~ 4 ± 1 and 
~ 2.5 ± 1 respectively. Both models yield a good match to the ESS BVR C redshift distributions to 21 — 22 mag and 
to the number-counts to 23 — 23.5 mas , which probe the galaxy distribution to redshifts z ~ 0.5 and z ~ 1.0 
respectively. 

The present analysis shows the usefulness of the joint use of the magnitude and redshift distributions for studying 
galaxy evolution. It also illustrates how Gaussian+Schechter composite luminosity functions provide more robust 
constraints on the evolution rate than pure Schechter luminosity functions, thus emphasizing the importance of 
performing realistic parameterizations of the luminosity functions for studying galaxy evolution. 

The detected density evolution indicates that mergers could play a significant role in the evolution of late-type 
Spiral and dwarf Irregular galaxies. However, the ESO-Sculptor density increase with redshift could also be 
caused by a ~ i ma s brightening of the Sc+Sd+dl galaxies at z ~ 0.5 and a ~ 1.5 — 2.0 mag brightening at 
z ~ 1, which is compatible with the expected passive brightening of Sc galaxies at these redshifts. Distinguishing 
between luminosity and density evolution is a major difficulty as these produce the same effect on the redshift 
and magnitude distributions. The detected evolution rate of the ESO-Sculptor Sc+Sd+dl galaxies is nevertheless 
among the range of measured values from the other existing analyses, whether they provide evidence for density 
or luminosity evolution. 

Key words, galaxies: luminosity function, mass function - galaxies: evolution - galaxies: distances and redshifts - 
galaxies: spiral - galaxies: irregular - galaxies: dwarf 
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1. Introduction 

Since the availability of the deep optical numbers 
counts, the excess at faint magnitudes has provided 
the majo r evidence fo r galaxy evolution at increasing 
redshifts l|Tvsonl Il988t iLillv et all Il99lt iMetcalfe et all 
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119951) . Using mo dels of the spectro-photometric evolu- 
tion of galaxies llGuiderdoni fc Rocca-Vohnerangel ll990: 
iBruzual fc Chariot! I199.C either passive luminosity 
evolution or more complex effects have been sug- 
gested to explain the faint number-count excess 
llGuiderdoni fc Rocca-Volmerangell99lUBroadhurst et alJ 
Il992t iMetcalfe et al.lll995[) . Althoughthe excess objects 
where initially envisioned as bright early-type galaxies 
at high redshift, the lack of a corr esponding h igh red- 
shift tail in the redshift distribution l)Lilrvl ll993) consoli- 
dated the interpretation in terms of evolution of later type 
alaxies, namely Spiral an d/or Irregular/Peculiar galaxies 



IjCampos fc Sha nks 1997) 



Here, we report on yet another evidence for evolu- 
tion of the late-type galaxies, derived from the ESO- 
Sculptor Survey (ESS hereafter). The ESS provides a 
nearly complete redshift survey of galaxies at z < 0.5 
over a contiguous area of the sky. A reliable description 
of galaxy evolution require proper identification of the 
evolving galaxy populations and detailed knowledge of 
their luminosity functions. In this context, the ESS sam- 
ple has the advantage to be split into 3 galaxy classes 
whic h are based on a template- free spectral classifica- 
tion ijGalaz fc de Lapparentlll998|) . and which are domi- 
nated by the giant morphologic al types E+SO+Sa, Sb+Sc, 
and 



Sc+Sd+Sm respectively l|de Lapparent et alJ 
Paper I hereafter). In Paper I, we have performed a de- 
tailed measurement of the shape of the luminosity func- 
tions (LF hereafter) for the 3 ESS spectral classes. The 
spectral-type LFs show marked differences among the 
classes, which are common to the B, V, R c bands, and 
thus indicate that they measure physical properties of the 
underlying galaxy populations. 

The analysis of the ESS LFs in Pap er I also pro- 
vides a revival of the view advocated by iBinggeli et alJ 
l)l988l) : a galaxy LF is the weighted sum of the intrinsic 
LFs for each morphological type contained in the con- 
sidered galaxy sample; in this picture, differences in LFs 
mark variations in the galaxy mix rather than variations 



1984 IBinggeli et alJ Il99fl 1 


Ferguson & Sandage 


1991: 


Trentham & Hodgkinl 2002: 


iTrentham & Tullvl 12002). 



Local measures show that giant galaxies (Elliptical, 
Lenticular, and Spiral) have Gaussian LFs, which are 
thus bounded at both bright and faint magnitudes, 
with the Elliptical LF skewed towards faint m agnitudes 
I Sandage et"aL1ll985t I.Terien fc Tammannl fl997^ . n con- 
trast, the LF for dwarf Spheroidal galaxies may be ever 
increasin g at faint magnitude s to the limit of the existing 
surveys ilSandage et alJll985t iFerguson fc SandageHl991t 



.erguson 

Jerien et alJ 1200ft iFlint et all l2001hlat IConselice et all 




, whereas the LF for dwarf Irreg ular galaxies is flat- 



ter ijPritchet fc van den Berghl Il999t) and may even be 
bounded at faint magnitudes IIFe rguson fc SandageH l989: 
Ijerien fc TammanrJll997t I Jerien et alJl2000|) . In Paper I, 
by fitting the ESS spectral-type LFs with composite func- 
tions based on the Gaussian and Schechter LFs measured 
for each morphological type in local galaxy groups and 



clusters ijSandage et al.l fl985t I Jerien fc Tammannl IT997fl . 
we confirm the morphological content in giant galaxies of 
the ESS classes, and we detect an additional contribution 
from dwarf Spheroidal (dE) and dwarf Irregular galaxies 
(dl) in the intermediate-type and late-type classes respec- 
tively. We then suggest that by providing a good match 
to the ESS spectral-type LFs, the local intrinsic LFs may 
extend to z ~ 0.5 with only small variations. 

In the following, we report on the measurement of the 
amplitude of the LFs for the 3 ESS spectral-type LFs, and 
on the detection and measurement of redshift evolution for 
the late-type galaxies. Sect. El lists the main characteris- 
tics of the ESS spectroscopic survey. Sect. El recalls the 
definition of the ESS spectral classes and the technique 
for deriving the corresponding K-corrections and absolute 
magnitudes. Sect. ^ shows the measured composite fits of 
the ESS LFs for the 3 spectral classes. In Sect. El we de- 
scribe the various techniques for measuring the amplitude 
of the LF (Sect.OOJ an d the associated errors fSect. l5~^|l : 
we then apply these techniques to the ESS and show the 
detected evolution in the late- type galaxies (Sect.EU- In 
Sect. El we use the ESS magnitude number-counts to de- 
rive improved estimate of the late-type galaxy evolution 
rate in the B, V, and R c bands. We then examine in Sect.0 
the redshift distributions for the 3 spectral classes in the 
3 filters, and we verify that the measured evolution rates 
for the late-type galaxies match the ESS expected redshift 
distributions. Then, in Sect. El we compare the detected 
evolution in the ESS LF with those derived from other ex- 
isting redshift surveys which detect either number density 
evolution fSect. l8~T|l or luminosity evolution fSect. 
Finally, Sect. El summarizes the results, discusses them in 
view of the other analyses which detect evolution in the 
late-type galaxies, and raises some of the prospects. 

2. The ESS spectroscopic survey 

The ESO-Sculptor Survey (ESS hereafter) provides a com- 
plete photometric and spectroscopic survey of galaxies in 

a region centered at ~ h 22 m (R.A.) 30°06' (DEC.), 

near the Southern Galactic Pole. The photometric sur- 
vey provides standard magnitudes B, V , and R c in the 
Johnson-Cousins system, for nearly 13000 galaxies to V ~ 
24 over a contiguous rectan gular area of ~ 0.37 deg 2 
[1.53°(R.A.) x 0.24°(DEC.)] llArnouts et alJll997ft . The 
uncertainties in the apparent magnitudes are 0.05 mag in 
the B V and R c bands for Rc ^ 21.0 l(Arnouts et alJ 
Il997l) . Mult i-slit spectroscopy of t he ~ 600 galaxies with 
R c < 20.5 l|Bellanger et al.lll995|) have provided a 92% 
complete redshift survey over a contiguous sub-area of 
- 0.25 deg 2 [1.02°(R.A.) x 0.24°(DEC.)]. Additional red- 
shifts for ~ 250 galaxies with 20.5 < R c < 21.5 were 
also measured in the same sub-area, leading to a 52% red- 
shift completeness to R c < 21.5 (see Paper I for details). 
We also consider here the V < 21.0 and B < 22.0 red- 
shifts samples, which correspond to the combination of 
the R c < 20.5 "nominal" limit with the typical colors of 
galaxies at that limit: B — R c ~ 1.5 and V — R c — 0.5 
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ijArnouts et al.lll997|) '). The redshift completeness for the 
V < 21.0 and B < 22.0 samples is 91% and 86% respec- 
tively. 

3. Spectral classification and K-corrections 

Estimates of morphological types are not available for the 
ESS redshift survey. In Paper I, we estimate the ESS 
intrinsic LFs based on a spectral classification. Using a 
Principal Component Analysis, we have derived an ob- 
jective spectral sequence, which is parameterized continu- 
ously using 2 parameters, describing respectively the rela- 
tive fractions of old to young stellar populations (param- 
eter denoted here Ts), and the relative st rength of the 
emission lines (|Galaz fc de Lapparentlll998(l . 

The ESS spectral sequence is separated into 3 classes, 
denoted "early- type" , "intermediate-type", and "late- 
type", which correspond to T s < -5°, -5° < T s < 3°, 
and Ts > 3° respectively. These values separate the ESS 
spectroscopic R c , V, and B samples into sub-samples with 
as least 100 galaxies (see Table QJ. Given the moderate 
number of objects in the ESS spectroscopic sample, these 
3 classes thus provide a satisfying compromise between 
resolution in spectral-type and signal-to-noise in the cor- 
responding LFs. 

Projection of the Kennicutt (1992) spectra onto the 
ESS spectral sequence shows that a tight correspondence 
wit h the Hubble morphological s equence of normal galax- 
ies l|Galaz fc de Lapparendll998t Paper I); the ESS early- 
type class contains predominantly E, SO and Sa galax- 
ies, the intermediate-type class, Sb and Sc galaxies, and 
the late-type class, Sc, Sd and Sm/Im galaxies. We show 
that these spectral classes allow us to detect the respec- 
tive contributions to the LF from the Elliptical, Lenticular 
and Spiral galaxies, and from the dwarf Spheroidal and 
Irregular galaxies. 

In Paper I, we estimate the K-correction for the ESS 
by projecting templates extracted from the PEGASE 1 
spectrophotometric model of galaxy evolution onto the 
ESS spectral sequence. A fine mesh of model spectra with 
varying redshift and spectral-type are generated, and K- 
corrections in the B, V, and R c filters are calculated for 
each of them; because the model spectra extend from 2000 
A to 10000 A, we can derive K-corrections in the 3 bands 
up to z ~ 0.6, the effective depth of the ESS redshift sur- 
vey. A 2-D polynomial fit to the resulting surface in each 
filter provides analytical formulas for the K-corrections 
as a function of redshift and spectral type, which are 
used for the ESS galaxies; these are plotted in Fig. Elbe- 
low, where they are c ompared with the K-corrections by 
IColeman et alJ l)l980t see Paper I for further details and 
comparisons) . For each galaxy, the absolute magnitude M 
is then derived from the apparent magnitude m, the spec- 
tral type Ts and the redshift z using 

M(m, z,T s ) = m-5 \ogd h (z) - K(z, T s ) - 25, (1) 

1 Projet d'Etude des GAlax ies par Synthese Evolutive 
iFioc fc Rocca-Volmerangelll99 , it) . 



wher e dj J (z) is the luminosity distance in Mpc l|Weinberel 
Il976h : 

d h (z) = (i + z) -i- / [n ni (i + v ) 3 + n A ]-* dv (2) 

H o Jo 

Throughout the present analysis, absolute magnitudes are 
calculated with a Hubble constant at present epoch writ- 
ten as Ho = lOOh km s" 1 Mpc -1 . We also assume a flat 
Universe with values of the dimcnsionlcss matter den- 
sity and cosmological constant assign ed to fl m = 0-3, 
ri\ = 0.7 resp., as cu r rently favored llRiess et al . 1998; 
Perlmutter et alJ Il999t IPhillips et all l200lt iTonrv et alJ 
2003); when specified, Sl m = 1.0 and 57 a = 0.0 are also 
considered (see Sects, 0] and 

In Paper I, we report on all sources of random and sys- 
tematic errors which affect the spectral classification, the 
K-corrections, and the absolute magnitudes. By compar- 
ison of the 228 pairs of independent spectra, we measure 
"external" errors in the K-corrections from 0.07 to 0.21 
mag, and resulting uncertainties in the absolute magni- 
tudes from 0.09 to 0.24 mag (with larger errors in bluer 
bands for both the K-corrections and the absolute mag- 
nitudes). From the 228 pairs of spectra, we also mea- 
sure an "external" r.m.s. uncertainty in the redshifts of 
(j ~ 0.00055, which causes negligeable uncertainty in the 
absolute magnitudes compared to the other sources of er- 
ror. 

4. The shape of the ESS luminosity functions 

The shape of the LFs for the 3 ESS spectral classes are 
derived in Paper I, using both the non-parametric step- 
wise maximum likeliho od method (SWML ) developed by 
lEfstathiou et alJ l|l988|) . and the method of lSandaee et alJ 
ill979l denote STY) which assumes a specific parametric 
form for the LF. Although pure Schechter (1976) functions 
provide acceptable STY fits to the 3 ESS spectral-classes, 
we show that as good or better STY fits of the R c LFs are 
obtained using composite functions based on the intrinsic 
LFs per morphological type measured in local groups an d 
clusters ijjerien fc Tammannlll997t ISandage et al.lll985[) . 

For the ESS intermediate-type and late-type R c LFs, 
we fit the sum of a Gaussian component representing the 
giant galaxies (Sb+Sc, and Sc+Sd resp.), and a Schechter 
component representing the dwarf galaxies (dSph and dl 
resp.). The Gaussian component is parameterized as 

<f>(M) dM = o e- (M °- M)2/2s2 dM, (3) 

where Mq and S are the pea k and r.m.s. dispersion respec- 
tively The ISchechterl l|l976j) component is parameterized 
as 

0(L)dL = ^(£)V*d(il) (4) 

where <jf is the the amplitude, L* the characteristic lu- 
minosity, and a determines the behavior at faint lumi- 
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Table 1. Parameters of the Gaussian and Schechter components of the composite luminosity functions fitted to the 
ESO-Sculptor spectral- type luminosity functions, in the R c , V, and B filters. 



Sample 


Numb. 


Morphol. 


Gaussian component 


Schechter component 


4>o 




of gal. 


content 


M -51ogh E a E t 


M* — 5 log h 


a 


0.4 In 10 (f>* 


early-type galaxies 


R c < 21.5 


291 


E+SO+Sa 


-20.87 ±0.23 0.84 ±0.24 1.37 ±0.36 








V < 21.0 


156 


E+SO+Sa 


-20.37 0.84 1.37 








B < 22.0 


108 


E+SO+Sa 


-19.47 0.84 1.37 








intermediate-type galaxies 


R c < 21.5 


270 


Sb+Sc/dSph 


-20.27 ±0.21 0.91 ±0.18 


-19.28 ±0.37 


-1.53 ±0.33 


0.83 


V < 21.0 


169 


Sb+Sc/dSph 


-19.87 0.91 


-18.88 


-1.53 


0.83 


B < 22.0 


154 


Sb+Sc/dSph 


-19.27 0.91 


-18.28 


-1.53 


0.83 


late-type galaxies 


R c < 21.5 


309 


Sc+Sd/dl 


-19.16 ±0.29 0.97 ±0.13 


-18.12 ± 0.22 


-0.30 


0.1 


V < 21.0 


168 


Sc+Sd/dl 


-18.76 0.97 


-17.72 


-0.30 


0.1 


B < 22.0 


190 


Sc+Sd/dl 


-18.36 0.97 


-17.32 


-0.30 


0.1 



Notes: 

- For the 2-component luminosity functions fitted to the intermediate-type and late-type galaxies, the morphological content 
of the 2 components appear separated by a "/" . 

- The listed luminosity function parameters correspond to Hq = 100ft km s _1 Mpc" 1 , Q m = 0.3, and Qa = 0.7. 



nosities. Rewritten in terms of absolute magnitude, EqJ3| 
becomes: 

4>{M) dM = 0.4 In 10 <f>*e- x X a+1 dM 

with (5) 

X = L_ = iqQA{M"-M) 

where M* is the characteristic magnitude, and a + 1 the 
"faint-end slope" . 

For the ESS early-type R c LF, a two-wing Gaussian 
function is used (a Gaussian with two different dispersion 
wings at the bright and faint end), as it successfully re- 
flects the combination of a skewed LF towards faint magni- 
tudes for the Elliptical galaxies with a Gaussian LF for the 
Lenticular galaxies displaced towards brighter magnitudes 
and with a narrower dispersion. The two-wing Gaussian 
is parameterized as 

4>{M) dM = Qe -(Mo-M) 2 /25^ dM foj . M < Mo 

= fae-iMo-Mf/svl dM for M > Mq W 

where Mq is the peak magnitude, and S a and S& are the 
dispersion values for the 2 wings. The improvement of the 
composite LFs over pure Schechter LFs is most marked for 
the early-type LF, because of its bounded behavior at faint 
magnitudes which is better adjusted by a Gaussian than a 
Schechter function, and for the late-type LF which cannot 
be fitted by a single Schechter function at simultaneously 
faint and intermediate magnitudes. 

The ESS R c LFs are measured for both the R c < 20.5 
and R c < 21.5 samples in Paper I. Here we choose to use 
the LFs measured from the deeper sample because the 
faint-end of the LF is better defined than from the shal- 
lower sample. Composite fits to the LFs for the V < 21.0 



and B < 22.0 samples have not been performed in Paper I. 
Here, we however need these fits in order to use the con- 
straints on galaxy evolution provided by the V and B faint 
number counts. Instead of performing the STY composite 
fits in the V and B bands, which requires particular care 
because of the number of parameters involved, we prefer to 
use the simpler approach of converting the R c LF param- 
eters into the V and B bands. We showed in Paper I that 
the shift in magnitude of the V and B LFs with respect to 
the R c is close to the mean V— R c , B—R c resp. color of the 
galaxies in the considered spectral class. Moreover, the se- 
lection biases affecting the V and B samples (a deficiency 
in galaxies bluer than B — R c ~ 1.5 and V — i? c — 0.5 
resp., due to the selection of the spectroscopic sample in 
the R c band) affect the faint-end of the late- type LF when 
fitted by pure Schechter functions, making it flatter than 
in the R c band. Conversion of the composite R c LF into 
the B and V bands allow us to circumvent the problem of 
incompleteness in the B and V bands. 

Figure plots the LFs for the 3 galaxy types in 
the R c < 21.5, V < 21.0 and B < 22.0 samples. 
The points represent the SWML solutions derived with 
(n m ,n A )=(l-0,0-0) in Paper I. For the R c < 21.5 sample 
(left panels) , the curves show the composite fits derived in 
Paper I; note that for the late-type LF, we have adopted 
the intermediate slope a — —0.3 (between the values 
a = —0.8 and a = 0.39 measured from the R c < 20.5 and 
R c < 21.5 samples; a slope a ~ —0.3 is also measured for 
Sm/I m galaxies in the Virgo cluster l|jerien fc Tammannl 
ll997l) L The parameters for the V and B composite LFs 
(middle and right panels resp. of Fig. are then de- 
rived from those for the R c < 21.5 sample by applying 
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M(RJ - 5 log h M(V) - 5 log h M(B) - 5 log h 

Fig. 1. The ESO-Sculptor Gaussian+Schechter composite luminosity functions for the early-type (top panels), 
intermediate- type (middle panels), and late- type galaxies (bottom panels) in the 3 samples: R c < 21.5 (left pan- 
els), V < 21.0 (central panels) and B < 22.0 (right panels). Symbols indicate the SWML solution, and lines the 
composite fits obtained; all were obtained with Ho = lOO/i km s _1 Mpc -1 , Q m = 1.0, and fl\ = 0.0. The shape 
parameters of the composite functions converted to Q m = 0.3 and = 0.7 are listed in Tabled 



the mean M(V) - M(R C ) and M (B) - M(R C ) colors for 
each spectral class (we use the R c < 20.5 sample for the 
color estimation, rather than the R c < 21.5 sample, as 
the completeness at magnitudes fainter R c = 20.5 is bi- 
ased in favor of red objects). The colors are those listed 
in Table 4 of Paper I: M{V) - M{R C ) = 0.5,0.4,0.4 and 
M(V)-M(R C ) = 1.4, 1.0, 0.8 for early-type, intermediate- 
type, and late-type galaxies respectively. 



Note that the SWML points in Fig. ^ account for the 
incomplete ness per apparent magnitude interval, as de- 
scribed bv IZucca et"aLl l|l994|) . For the SWML points, a 
bin size of AM = 0.48 mag is used in all filters (smaller 
or larger bin sizes within a factor 2 yield similar curves). 
As the amplitudes of the composite fits and the SWML 
solutions in the Fig. ^ are so far undetermined (they are 
measured later on in Sects. [S] and EJ, we adopt the follow- 
ing: we use the same normalization of the SWML curves 
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as used in Paper I (see Table 3), using the amplitude </>* 
measured from the pure Schechter fits; then for each sam- 
ple, the composite is adjusted by least-square fit to the 
SWML points (with the ratio </>o/0.4 In 10 0* between the 
Gaussian and Schechter component kept fixed to the val- 
ues in Table Pi . 

Figure ^ shows that the composite LFs provide good 
adjustment to the SWML solutions for each of the 3 spec- 
tral classes in each filter. In particular, the simple color 
shift used to define the composite LFs in the V and B 
bands provides good adjustments to the SWML points in 
both bands, despite the color biases affecting the redshift 
completeness of these samples. Note that the composite 
spectral- type LFs derived from the R c < 21.5 sample pro- 
vide satisfying adjustment to the SWML points for the 
R c < 20.5 spectral-type samples. 

In Paper I, the LFs are measured for cosmologi- 
cal parameters f2 TO = 1.0 and f^A = 0.0. Because the 
faint number counts which we use below to constrain 
the ESS evolution rate are sensitive to the cosmologi- 
cal parameters, we have converted these values to O m = 
0.3 and fl\ = . 7, the currently favor e d parameters 



llRiess et al.lll998t IPerlmutter et all Il999t IPhillips et all 
200^rT^nrvet"ld1 12003|) . Again, rather than re-running 



the com posite fits, we apply the empirical corrections de- 
rived bv lde LapparentJ (|2T 



as follows. When changing 
from (O ro ,O A )=(Q.3,0.7) to (fl m , fi A )=(l-0,0.0), the vari- 
ation in absolute magnitude due to the change in luminos- 
ity distance is AM ~ 0.3 mag at z ~ 0.3, the peak redshift 
of the ESS (see Figs. 113115(1 . T his empirica l corre ction 
is confirmed by the results from iFried et al.l l|200ll) and 
iBlanton et al.1 ll200ll). who calcu late galaxy LFs in both 
cosmologies. Ide Lapparentl l|2003j) also apply a correction 
to the Schechter parameter a, with 8a ~ AM/3, due to 
the strong correlation between the M* and a parame- 
ters in a Schechter parameterization. Here we neglect this 
correction, which would amount to 5a ~ 0.1 for a pure 
Schechter parameterization; this value is comparable or 
smaller than the l-er uncertainty in the faint-end slope 
a of the Schechter component for the ESS intermediate- 
type and late-type LFs, and than the 1-er uncertainty in 
the dispersion £ a and of the 2-wing Gaussian fitted to 
the ESS early- type LF (see Table [IJ. 

Tablc^lists the resulting LF shape parameters for the 
Gaussian and Schechter components (Mo, E a , £{,, M*, a; 
for a symmetric Gaussian, £ is listed in the £ Q column): 
conversion to (Q m ,fiA)=(0.3,0.7) is obtained by shifting 
all values of M and M* by -0.3 mag . The last column of 
TableQJlists the ratio of amplitude cj>o /OA In 10 <fi* between 
the Gaussian and Schechter component derived from the 
composite fits to the R c < 21.5 sample in Paper I. We 
adopt the same values of this ratio for the V and B LFs. 

5. The amplitude of the ESS luminosity functions 

The amplitude of the LF is proportional to the mean den- 
sity of galaxies brighter than some absolute magnitude 
threshold. It is therefore a useful indicator of the large- 
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Fig. 2. Absolute R c magnitudes as a function of redshift z 
for the ESS galaxies with R c < 21.5, using Ho — lOOh km 
s _1 Mpc -1 , Cl m — 0.3, and FIa — 0.7. This graph shows 
how the apparent limiting magnitude biases the range of 
absolute magnitudes detected at increasing redshift, and 
how the varying K-corrections per spectral type affect the 
faintest absolute magnitude reached at a given redshift. 



scale variations in the luminous component of the matter 
density in the Universe, and its possible evolution with 
redshift. 

Several other redshift surveys probing the same red- 
shift range as the ESS survey (0.1 ^ z ^ 0.5) have de- 
tected signs of evolution in the intrinsic galaxy LFs (see 
Sect.[H]below). The major difficulty in measuring evolution 
of the LF with redshift originates from the limit in appar- 
ent magnitude which affects most redshift surveys. The 
flux limit results in the detection of galaxies in an absolute 
magnitude range which narrows with increasing redshift 
to the brightest galaxies. Figure [21 shows how the absolute 
magnitudes of the faintest detected ESS galaxies is a func- 
tion of redshift and spectral- type (via the K-correction) : 
the limiting curve at faint magnitudes is defined by replac- 
ing in Eq.r^the apparent magnitude m with the R c — 21.5 
magnitude limit, and the K-correction K(Tg,z) with the 
smallest value over the considered spectral class at each 
z. Faint galaxies (—18 + 5 log h ^ MrR c ) ^ —15 + 5 log h) 
are exclusively detected at z ^ 0.2 in the ESS, whereas 
only bright early-type and intermediate- type galaxies with 
M(R C ) < —21.0 + 5 log ft,, and bright late-type galaxies 
with M(R C ) < -20.0 + 5 log ft are detected at z > 0.5. 
In the full ESS redshift range 0.1 ^5 z ^ 0.5, only galax- 
ies in the magnitude interval —22.0 + 51og/i ^ Mr 
—21.0 + 5 log ft. can be observed, which is clearly too nar- 
row for deriving any constraint on the evolution in the 
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Fig. 3. Same as Fig. for V < 21.0. 



Fig. 4. Same as Fig.|2jfor B < 22.0. 



shape of the LF. Figures [3] and 0] which show the abso- 
lute V and B magnitudes versus redshift for the 3 spectral 
classes in the V < 21.0 and B < 22.0 samples resp., dis- 
play similar effects. IPriverl l|200lll showed that these biases 
strongly affect the usual tests for galaxy evolution based 
on the shape of the LF. 

The other limitation for detecting redshift evolution 
in the ESS is the limited statistics: separation of each of 
the 3 spectral classes into even as few as 2 redshifts inter- 
vals would yield large uncertainties in the measured shape 
of the LFs, which would make insignificant any reason- 
able difference between the high and low redshift LFs. For 
the ESS, we can only examine whether the "general" LF, 
i.e. the LF summed over all ESS spectral types, evolves 
with redshift. Here we consider the LF at R c < 20.5, 
as varying incompleteness at fainter magnitudes may act 
as evolution. For (Q m ,fiA.) = (1-0,0.0), the ESS general 
LF for R c < 20.5 can be fitted by a Schechter function 
with M* = -20.95 ±0.12 + 5 log h and a = -1.15 ±0.09 
(STY fit). The corresponding LF in the redshift interval 
z < 0.3 at R c < 20.5 has M* = -21.17±0.21 + 51og/i and 
a = -1.18±0.11; at z > 0.3, M* = -21.17±0.28 + 5 logft 
and a = —1.63 ± 0.31. The values of M* are identi- 
cal in the 2 redshift ranges, and the slope a differs by 
1.4-sigma; using (O m , fiA)=(0.3,0.7) yields similar conclu- 
sions. Given the difficulty in determining the parameter 
a in a Schechter fit (Paper I), these results are consistent 
with the hypothesis of null evolution in the shape of the 
ESS Rc "general" LF in the redshift range 0.1 <, z < 0.5. 
This is however no proof that the individual spectral-type 
LFs do not vary in shape: one could imagine redshift vari- 
ations in the shape of the intrinsic LFs which would con- 
spire to combine into a constant "general" LF. 



IPriverl l|200lh did suggest that reliable constraints on 
evolution of the shape of the galaxy LF at z ^ 1.25 
may only be derived from sample s as deep as the Hubble 
Deep Field (Wi lliams et~ai1ll 996). which reaches z ~ 3.5 
( Sawicki et al.ll997tl : lDriverl l|200l(l recommends that such 
analyses be rather based on the bi-variate brightness dis- 
tribution for galaxies (the function which describes the 
galaxy bi-variate distribution in absolute magnitude and 
mean surface brightness). To circumvent the difficulty in 
measuring evolution in the shape of the LF, we assume in 
the following that the shape of the intrinsic LFs for the 
ESS early-type, intermediate-type, and late-type galaxies 
(as listed in Table ^1 is not evolving with redshift. Any 
possible evolution will then be detected as evolution in 
the amplitude of the LF. 

5.1. Measuring the amplitude of the luminosity 
function 

Once the shape of the LF is determined, its amplitude 
can be determined in a second stage. We separate the 
amplitude $ from the "shape" component <p{M) dM: 

<p(M) dM = $ tp(M) dM. (7) 

For a Gaussian LF, $ = <po, and for a pure Schechter LF, 
$ = 0.4 In 10 4>* (see Eqs. 01 and 15- In a survey where 
the detected brightest and faintest absolute magnitudes 
are M\ and Mi resp., one can calculate the mean density 
of galaxies with Mi < M < M 2 , denoted n (Mi < M < 
M%), which is related to <& by 

n (Mi < M < M 2 ) , , 

$ = — !: — _ = — -L- r«l 

C„(M)dM ' [) 
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In a magnitude-limited survey, the absolute magni- 
tudes Mbright(z) and Mf aint (z) of the detected brightest 
and faintest galaxy resp. vary with redshift (see Figs. |2J 
E}- Calculation of the mean density n (Mi < M < M 2 ) 
therefore requires to correct the observed number of galax- 
ies to the expected number if the survey was limited to the 
constant absolute magnitude interval Mi < M < M 2 . The 
correction is obtained by multiplying the observed num- 
ber of galaxies at redshift z by the inverse of the selection 
function S(z) defined as 



S(z) 



,min(M f . in t(*) Aft) (M) dM 



AM 



(9) 



For the ESS, we take in all filters M\ — —30 + 5 log h, and 
M 2 = —16 + 5 log h (see Figs. 12141) . S(z) then measures 
the fraction of galaxies with Mi < M < M 2 at redshift z 
which are inc luded in the survey. 

Following iDavis fc Huchral lfl982ft . we define 3 estima- 
tors for the mean density n (Mi < M < M 2 ) which 
are unbiased by the apparent magnitude limit of the 
survey (in the following, although we omit to mention 
Mi < M < M 2 , all quoted densities refer to that interval). 
If N(z) dz is the observed number of galaxies in a shell 
dz at redshift z, N(z) dz/S(z) is the expected number 
of galaxies with Mi < M < M?. A first est i mator of the 
mean density is defined bv lDavis fc Huchral l|l982ft as 



fll (-^max) 



1 



N(z) 



dz, 



(10) 



V( 

where z m j n and z max are arbitrary choices of the small- 
est and largest redshift over which the integrals are per- 
formed, and V(z m - m , z max ) is the total volume of the survey 
between these limits: 



V(z 



mm ; v max J 



— — dz 

dz 



(11) 



(dV is the comoving volume element at redshift z, see 
Eq.[23J). Here we use z m [ n = 0.1 (see Sect. I5.3f) . Note that 
because 1/S(z) rises sharply with reds hift, the n-\ esti- 
mator heavily weights distant structures. IDavis fc Huchral 
l)l982j) also showed that ni is close to the minimum vari- 
an ce estimator of th e mean density. 

IDavis fc Huchral l)l982j) define another estimator by 
equating the observed number of galaxies with the ex- 
pected number in a homogeneous universe: 



^3 (^max) 



s:::s(z)%dz 



(12) 



This estimator is much more stable than n\, as all ob- 
served galaxies are weighted equally. It however heavily 
weights galaxies near t he peak of the redshift distribution 

foavis fc Huchralll982l). 

Finally. iDavisfe" Huchral l|l982() define a third and in- 
termediate estimator by averaging the expected density 
across radial shells: 



l2(2max) 



N(z) 



dz 



(13) 



The interest of the n 2 estimator is that its differential 
value 



dm 



N(z) 



dz 



(14) 



can also be calculated as a function of redshift, and al- 
lows one to examine the variations of the local galaxy 
density with redshift. Note that the 3 estimators ni, n 2 , 
and 713 can be calculated to a varying depth z max and are 
thus defined as functions of z max . In a fair sample of the 
galaxy distribution, that is a sample which is significantly 
larger than the largest density fluctuations, all 3 estima- 
tors should converge to a common value as z max reaches 
the sample depth. 

A fou rth estimator, which w e test here, is that pro- 
posed bv lEfstathiou et all l|l988j) . This estimator is simi- 
lar to ni but the shell dz is taken to be infinitely small so 
that galaxies can be counted one by one, and the integral 
can be re- written as a sum over the galaxies in the sample: 



"-EEp(Zmax) — 



1 



1 



V(z 



mini -° in ax. 



(15) 



where Zj is the redshift of each galaxy, JV ga i galaxies are 
observed in the redshift range z m j n — z max , and the volume 
^(zmin, ^max) is defined in Eq.[TTl 

The variations with z max of the estimators $i, $2, $3, 
and $EEP, based on ni(z max ), n 2 (z max ), n 3 (z max ), and 
«EEp(zmax) respectively, can then be defined using Eq.|H] 
In practise, the integrals in the m, n 2 and 713 estima- 
tors are calculated as discrete sums over a finite redshift 
bin Az, and the selection function correction is approxi- 
mated as S(z c ) where z c is the central redshift of the bin. 
However, because the second derivative of the selection 
function is positive, this yields an underestimate of the 
expected number of galaxies. The resulting systematic er- 
ror in $ is ^ 1%. It is reduced to ~ 0.1% by replacing 
S(z c ) by its average value over the bin, which we adopt 
for the nl, n 2 and estimators. The tt-eep estimator is 
a priori unbiased by this effect, as galaxies are considered 
one by one, and S(z) is calculated at the redshift of each 
galaxy. Our tests with mock ESS catalogues (described 
in the next Sect.) however show that the 71eep estimator 
tends to over-estimate the true $ by ~ 1% for uniform 
distributions with ~ 200 points; this bias disappears for 
distributions with more than ~ 2000 points. In the ESS 
spectral classes for which A^ ga i ~ 100 — 300 galaxies, a 
*~ 1% bias is small compared to the random uncertainties 
in the nl, n 2 and 723 estimators (see next Sect.). 

5.2. Estimation of errors 

We estimate the random errors in the amplitude of the 
LFs for the ESS by generating mock ESS distributions 
with - 240, ~ 2400 and ~ 24 000 points, and a Schechter 
LF defined by M* = —19.2 and a = —1.4 (using com- 
posite LFs as those listed in Table ^ would not change 
any of the reported results); the measured values of <fi* 
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are then in proportion of the total number of galaxies in 
each simulation. The mock distributions have no built-in 
clustering, and are thus uniform spatial distributions mod- 
ulated by the selection function (see Eq. . To neverthe- 
less take into account the effects of large-scale clustering, 
we also introduce in the simulations density fluctuations 
in redshift which mimic large-scale structure: fluctuations 
resembling those measured in the ESS (measured as the 
departure from the redshift distribution for a uniform dis- 
tribution) , and various transformations of these including 
the inverse of the ESS fluctuations. We neglect density 
fluctuations transverse to the line-of-sight, as they are ex- 
pected to have a ncgligcablc effect on the calculation of 
the mean density. 

For each type of simulation (choice of <$* and fluctua- 
tions in density), we generate 100 samples with different 
seeds for the random generators. The random errors in the 
various estimators of cf>* are then calculated directly from 
the variance among the 100 realizations of the same distri- 
bution, and arc functions of the redshift z max out to which 
Eq. 1 1011 131 are integrated. We measure that the variance 
in the estimator $3 of <f>* (based on 713) varies as 



g($3 

$3 



1.0 



mini ^max 



(16) 



where -ZV ga i(z m i n , z max ) is the number of galaxies in the 
survey which lie in the redshift interval [z m j n , z max ]- This 
is valid out to the most distant galaxy included in the 
survey. For the ni estimator of <p* , we measure 



1.5 



(17) 



This is also valid out to the redshift where N ga \ cor- 
responds ~ 95% of the total number of galaxies in 
the sample (z ~ 0.6 in the ESS). For the m and 
7jeep estimators, the relative variance in </>*(z zmax ) rise 
slowly from 1.0/ \/iVg a i(.z m i n , z max ) at low redshift to 
1.5/ ^/TVga^Zmin, z max ) when TVgai reaches ~ 85% of the to- 
tal number of galaxies in the sample (z ~ 0.4 in the ESS) 
and can grow to 2.0/y / -/V ga i(;z m i n , z max ) at larger distances 
(in the redshift interval 0.4 < z < 0.6 for the ESS). 

Davis & Huchral |l982) calculate that for the mini- 
mum variance estimate of the mean density (close to rii), 
the relative uncertainty due to galaxy clustering within 
the finite volume sample is 



a(n) 



min) *max 



V(z 



mm; ^max 



(18) 



where 47rJ3 is the volume integral of the 2-point galaxy 
correlation in the survey volume bounded in redshift by 
[z m i n , z max ] . By using th e correlation function £(rp,7r ) 
measured from the ESS bv lSlezak fc de Lapparent! l|2004|) . 
we calculate the values of J3(z m ; n , z max ), and we confirm 
that the relative uncertainties derived from Eq.^Jfor one 
spectral class are comparable to those given by Eq.ll6lfor 
Zmin = 0.1 and 0.1 < z max < 0.6, with Eq. 1181 becoming 
slightly larger than Eq. ^]at z & 0.3. 



As pointed bv lLin et all l|l996|) . Eq.UJand the errors 
calculated from the mock distributions underestimate the 
error in <jf , as they do not take into account the uncer- 
tainty in the shape of the LF: indeed, in our mock distri- 
butions, the LF is not re- measured, as we assume that the 
true value is known. In that case (when the input values 
of M* and a are used to calculate </>*), all 4 density es- 
timators recover the true value to 0.1 — 1% (see Eqs. [TBI 
and I17[). The systematic un derestimation of <fi* by 20% 
calculated bv lWillmerl l(l997h for LFs derived by the STY 
estimator, could be a measure of this effect, if one assumes 
that <fi* was measured using the biased values of M* and 
a (but the author unfortunately does not specify). Here, 
however, we do not need to evaluate this additional source 
of uncertainty because the large amplitude of the density 
fluctuations on scales of 100 hr 1 Mpc which are present 
in the ESS survey cause systematic variations in the deter- 
mination of both <j>* and <fio (see Sect. 15.31 which dominate 
all the other sources of error discussed here. 



5.3. Evolution in the ESS mean density 

Figure [5] shows for each of the 3 ESS spectral classes the 
variations with redshift z max of the estimates <!>i, $2, $3, 
and $eep of the amplitude 0o for the Gaussian component 
of the composite LFs (see Tableland Eqs. 01 and EJ); these 
estimates are based on the ni(z max ), n 2 (z m ax), n 3 (z mlix ), 
lEEp(zmax) estimators of the mean density resp. (see 
Eas. fTTMB) for the ESS R c < 21.5 sample. All quantities 
plotted in Fig. (and the other figures in this section) use 
Zmin = 0.1, in order to avoid biasing of the density esti- 
mates, as the small volume probed at z < 0.1 suffers from 
under-sampling (see Figs. l2l4|) . The estimators n\, ri2, n$ 
are calculated in steps of Az = 0.02. The random errors 
due to the finite sample size for the n\, 713, and 71eep es- 
timators are evaluated using Eq.^l f° r the 712 estimator, 
Eq. El is used. 

We emphasize that under the assumption of a non- 
evolving shape of the intrinsic LFs <p(M) for the various 
galaxy types (see Sect.EJ, examining the variations in </> 
with redshift is equivalent to examining the variations in 
the ESS mean density with redshift (see Eq.|BJ). Moreover, 
in the following, we only show and discuss the variations 
in the amplitude </>o of the Gaussian component of each 
spectral-type LF. Nevertheless, because the variations in 
the amplitude <fr* of the Schechter components for the 
intermediate-type and late-type LFs are simply propor- 
tional to those in </>o using the ratio 4>q/0A In 10 </>*, whose 
values are listed in Tabled all comments on the variations 
of 4>o also apply to those in 0*. In order to refer to both 
</>o and 0*, we use in the following the generic amplitude 

We first compare the performances of the 4 estimators 
ni , TI2, n3 and riREP- Figur eElshows that the tieep estima- 
tor ijEfstathiou et al.lll988j) yields nearly indistinguishable 
results from the n\ estimator. For all galaxy types, the 113 
estimator converges at z ^ 0.5 to a value determined by 
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Fig. 5. The amplitude </>o of the Gaussian component of 
the composite luminosity functions measured for the ESO- 
Sculptor R c < 21.5 sample, using the estimators ri\ (red 
solid line), (green dotted line), 123 (blue dashed line), 
and rtEEP (black dot-dashed line) defined in Eqs. HI1I IT51 
as a function of redshift (denoted z max in the text). Each 
panel corresponds to a spectral-type. For clarity, only 
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Fig. 6. The amplitude </>o of the Gaussian component 
of the composite luminosity functions measured for the 
ESO-Sculptor R c < 21.5 sample, using the estimator n\ 
defined in Eq. 1101 as a function of redshift z mia ; early- 
type, intermediate-type, late-type galaxies correspond to 
the red dotted, green dashed, and blue solid lines respec- 
tively. For clarity, the random errors for the intermediate- 



, . . type galaxies are not shown; they are comparable to those 

the random errors lor the ni estimator, calculated using . . . . , . , 

_ — . . , , ,• 1 , .1 \ mi 1 r lor the early- type and late-type galaxies (see Eg. 1161). ihc 



Eq.^|are shown (see text for details). The error bars for 
the other curves have similar amplitude, except for ri2, for 
which they are 50% wider. 



linear fit of Eq. ^5] is also plotted (magenta solid straight 
line). 



the ratio of total number of galaxies in the considered sam- 
ple by the integral of <p(M) (see Eq. 0); the asymptotic 
77,3 estimator of 4>o is however dominated by the observed 
number of galaxies at z < 0.3 (see Sect. 15 .l|) . In contrast, 
as the ri\ estimator weights different parts of the survey 
proportionally to their volume, it allows one to trace the 
large-scale density fluctuations across the ESS redshift in- 
terval. As expected, the 712 estimator of 4>q (Eq. 113(1 yields 
intermediate values between the n\ and 713 estimators for 
the 3 galaxy types. In the following, we thus restrict the 
discussion to the extreme cases represented by the n\ and 
n 3 estimators. 

FigurelHlprovides a direct comparison of the amplitude 
<fio of the Gaussian component of the LFs based on the 
ni estimator (denoted <j>i hereafter) for the 3 ESS spec- 
tral types. The 3 curves show a relative depression in the 
range 0.33 ^ z ^ 0.42 followed by an excess at z ~ 0.45. 



Both features correspond to large-scale fluctuations in the 
redshift distribution. Note that for the 3 galaxy types, the 
fluctuations in the $1 about its mean value have a typical 
amplitude of 0.001, corresponding to a relative variation 
of ~ 30%. These systematic errors largely dominate the 
random errors in $1 indicated by the vertical error-bars in 
Fig. (see also Eas.llb1and ll7|) . The other marked effect in 
Figgis the steady increase of <&i for the late-type galaxies 
compared to the more stable behavior for the early-type 
and intermediate-type galaxies: a factor ~ 2.5 increase is 
measured from z ~ 0.15 to z ~ 0.55. We interpret this 
effect as an evolution in $1 for the late-type galaxies and 
provide further evidence below. An increase of <i>3 for the 
late-type galaxies is also visible in Fig. [SJ but the effect is 
smaller as $3 applies a lower weight to the distant struc- 
tures relative to those nearby. 
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Table 2. n\ and n% estimators of the amplitude 0o of the Gaussian component of the early-type, intermediate-type, 
and late-type luminosity functions for the ESO-Sculptor redshift survey, in the Johnson B, V and Cousins R c bands. 



Sample 


Early- type 




Intermediate- type 


Late- type 










$i(0.51) 


$3(0.51) 


$i(0.51) 


$3(0.51) 






$i(0.51) 


$3(0.51) 


R c < 21.5 


0.00333 


0.00339 


0.00326 


0.00312 


0.00194(1 +3.51(2 


-0.15)] 


0.00439 


0.00328 


Rc < 20.5 


0.00324 


0.00331 


0.00351 


0.00321 


0.00177(1 +4.45(2 


-0.15)] 


0.00462 


0.00319 


V < 21.0 


0.00322 


0.00335 


0.00379 


0.00314 


0.00141(1 +7.52(2 


-0.15)] 


0.00521 


0.00309 


B < 22.0 


0.00286 


0.00363 


0.00357 


0.00325 


0.00135(1 +8.59(2 


-0.15)] 


0.00551 


0.00333 



Notes: 

- The listed values of the luminosity function amplitude <f>Q are in units of h 3 Mpc -3 mag -1 , and are obtained with Hq — 100/t 
km s _1 Mpc" 1 , n m = 0.3, and D. A = 0.7. 

- The corresponding amplitude 0* of the Schechter component for the intermediate-type and late-type galaxies can be derived 
using the values of the ratio 0o/O.41n 10 0* listed in Table Q 

- The $1(2) parameterization of the evolving amplitude <f>o for the late-type galaxies should be used with caution (especially 
in the V and B bands); more realistic parameterizations are obtained in the R c , V and B bands from the adjustment of the 
number counts (see Table 01 Sect.|SJ. 
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Fig. 7. The amplitude (f>o(z max ) of the Gaussian compo- 
nent of the composite luminosity functions measured from 
the ESO-Sculptor R c < 21.5 sample for the intermediate- 
type (green dashed line) and late-type galaxies (blue solid 
line), both normalized to 4>o(z maK ) for the early-type 
galaxies. All values of 4>o are derived using the n\ esti- 
mator defined in Eq. The linear fit to the relative late- 
type density given in Eq. I2J3 is also shown (magenta solid 
straight line). 



The variations of $1 with redshift in Fig. [B] thus re- 
sult from the combination of the density fluctuations pro- 
duced by the pattern of large-scale structure with the 
possible evolution in the amplitude of the LF. A sur- 
vey with a wider angular extent would be necessary to 
average out the effect of large-scale structures perpen- 
dicular to the line-of-sight. Although the different galaxy 
types have different clustering properties l|Lovedav et alJ 
Il995|) . which on small scales a re characterize d by the 
morphology-density relationship (Dressier 1980). they do 
trace the same patt ern of walls and voids at large scales 
illlucln-M ci a.lJ fl99f)1. There has been so far no detection 
of systematic variations in the proportions of the different 
galaxy types on scales larger than ~ 10 Mpc (and 
thus outside galaxy clusters). We therefore make the hy- 
pothesis that the proportions of the different galaxy types 
are constant at large scales, and we eliminate the fluctua- 
tions caused by the large-scale structure by normalizing $1 
for the intermediate and late-type galaxies by $1 for the 
early-type galaxies. The resulting relative variations of </>o 
in the R c < 21.5 sample are shown in Fig.0 The normal- 
ization erases most of the fluctuations in Fig. [BJ and only 
the deviations from a distribution having similar large- 
scale clustering as the early-type galaxies remain. Whereas 
the relative density for the intermediate-type galaxies re- 
mains within the narrow interval ~ 0.8 — 1.0 in the redshift 
interval 0.15 ^ z ^ 0.5, the relative density of late-type 
galaxies shows a linear increase by a factor of nearly 2 in 
this interval. 



Under the assumption that the shape (p(M) of the in- 
trinsic LFs for the various galaxy types are non-evolving 
with redshift (see Sect. [5J, the systematic increase with 
redshift in the $1 estimate of (j>o for the late-type galax- 
ies detected in Fig0 can be interpreted as evidence for 
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Fig. 8. The amplitude 0o( z max) of the Gaussian compo- 
nent of the composite luminosity functions measured from 
the ESO-Sculptor V < 21.0 sample (top panel) and the 
B < 22.0 sample (bottom panel) for the intermediate-type 
(green dashed line) and late-type galaxies (blue solid line) , 
both normalized to <^>o(-2max) for the early- type galaxies. 
All values of 4>o are derived using the n\ estimator defined 
in Eq. ^3 The linear fits to the relative late-type density 
and the corresponding parameterization are also shown in 
each panel (magenta solid straight line). 

evolution in the ESS. A linear regression to the relative 
late-type density in the interval 0.15 < z < 0.5 yields 



0cf tC 

.early 
y 



(z) = 0.582 [1 + 3.51 [z- 0.15)]; 



(19) 



this fit is over-plotted in Fig. [7| Note that the linear fit is 
performed in the restricted redshift interval 0.15 < z < 0.5 
because the survey volume is too small at z < 0.1 (see 
Figs. 121 II) and the redshift survey is too diluted at z > 
0.5 to provide a reliable estimate of 4>q. We now assume 
that <f>o for the early-type galaxies does not evolve with 
redshift at z < 0.5, and that the redshift evolution in 
Eq. E3 can be fully attributed to the late- type galaxies. 
The parameterization of the late-type density evolution 
therefore becomes 



-P0.15 
</>u ate 



= 3.51 

(0.15) = 0.00194 h 3 Mpc~ 



3 mag 1 



(21) 



The zero-point (/>jf tc (0.15) is measured by adjusting the 
linear fit to the measured values of the $i estimate of 
4>o at z = 0.51 for the late-type galaxies, which is listed 
in Table The linear evolution model of Eq. 123 is over- 
plotted in Fig. El with extrapolation to the redshift inter- 
vals 0.1 < z < 0.15 and 0.5 < z < 0.6. This model follows 
satisfactorily the increasing trend of the $i estimator of 4>q 
for the late-type galaxies (recall that the deviations from 
the model are mostly caused by the large-scale structures). 

We emphasize that the $i(z) estimator at redshift z 
is a cumulative measure over the redshift interval 0.1 — z, 
which may therefore underestimate the evolution rate for 
the late-type galaxies. However, the contribution at each 
redshift is proportional to the volume sampled; as the vol- 
ume increases with z 2 , this puts most of the weight on the 
structures at z, making $i a good approximation to an in- 
cremental measure of the density variations with redshift. 
We also confirm this result by measuring the estimate of 
<po using the incremental estimator of the mean density 
dri2, described in Sect. 15.11 (Eq. 114ft . The estimator d«2 
yields a similar rate of increase in <fio as the n\ estimator, 
but its large error bars prevent any reliable measure. 

Other authors have modeled the evo lution in the am- 
plitude $ of the LF as a power of 1 + z llLillv et al.lll996l 
iHevl et aT1ll997l) , as it converts to a po wer of cosmic tim e 
for (O m ,n A )=(1.0,0.0): 1 + z oc t~ 2 / 3 llCole et al.lll99l . 
Using the adopted values (f2 m ,f2A)=(0.3,0.7), adjustment 
of 0o atc / ( / ) o arly by a power-law function defined as 



6> ate (0) (1 + zF, 



(22) 



/late 
D 



(z) = ^ ate (0.15) [1 + Po.15 (z - 0.15)] 



(20) 



over the redshift interval 0.15 < z < 0.5, deviates from 
the linear model by at most ~ 0.0001 h 3 Mpc -3 mag -1 , 
significantly smaller than the 1-a random errors in the 
measurement of 0{f tc (see Fig. Therefore, in the lim- 
ited redshift range of the ESS, both the linear and power- 
law models provide good descriptions of the evolution of 
the relative density of late-type galaxies at R c < 21.5. In 
Sect. El we show that when extrapolated to z ~ 1, both 
the linear and power-law models can also adjust the ESS 
faint BVR C number-counts. 

For numerical comparison of the $i and $3 es- 
timators of (fto, we list in Table El their values at 
z max = 0.51 for each spectral class of the R c < 21.5, 
R c < 20.5, V < 21.0 and B < 22.0 samples. For 
the early-type and intermediate-type galaxies, the val- 
ues of $3(0.51) show a systematic difference of 0.00006 — 
0.0008 h 3 Mpc -3 mag -1 with the corresponding values of 
<i>i(0.51), due to the different weighting of the large-scale 
structure by the 2 estimators (see Sect. 15. l|) . In contrast, 
for the late-type galaxies, $3(0.51) is systematically lower 
than $i(0.51) by ~ 0.001 h 3 Mpc~ 3 mag -1 in R c to 
~ 0.002 h 3 Mpc~ 3 mag -1 in V and B. Ignoring the evo- 
lution in the late-type density and using the $3 estima- 
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tor of 0o rather than $i would then underestimate 0o at 
z ~ 0.5 by ~ 40 - 50% in i? c , and by - 70% in V and B. 

For the late-type galaxies, we also list in Table [21 
the linear parameterization of &i(z) (Eqs. 1201 and l?T|) . 
Comparison of the "zero-point" of the linear parameteri- 
zation <&i (0.15) with the listed value of <&i (0.51) illustrates 
the increase in <I>i from z ~ 0.15 to z ~ 0.51. Note that 
Table [21 also lists the parameters derived from the shal- 
lower but more complete R c < 20.5 sample (using the 
same LFs as for the R c < 21.5 sample, see Table QJ. The 
evolution rate P0.15 = 4.45 is close to that for the deeper 
R c sample, and the values of $1 and $3 differ from those 
from the R c < 21.5 sample by less than 1-a (the uncer- 
tainties in <!>! and $3 for all considered samples in Table[2] 
are in the interval ~ 0.0002 - 0.0003 for z max > 0.3). 

We apply a similar analysis to the ESS V < 21.0 and 
B < 22.0 samples. The measured and modeled relative $1 
estimators of 4>q are plotted in Fig. |H1 the relative <E>i for 
the intermediate-type galaxies are also plotted in Fig. 0D 
As in the R c filter, the fits are performed in the restricted 
redshift interval 0.15 < z < 0.5. In both the V and B 
bands, we find an increase of the late- type galaxy density 
which confirms the reality of the effect detected in the R c 
band. We however measure a factor of 2 higher evolution 
rate P0.15 in V and B than in R c (see Table EJ. This may 
be due to the fact that when going to bluer bands, the 
late-type galaxies are favored. The higher evolution rate 
in the V and B bands may also reflect the color-dependent 
selection effects affecting the redshift samples in the V and 
B bands. Moreover, the variations in fb.15 from band to 
band are symptomatic of the limited constraints provided 
by the ESS redshift distributions. In the following, we show 
that measuring the evolutions rates from the ESS faint 
magnitude number counts yields better agreement among 
the 3 filters. 



6. Modeling the ESS number-counts 

Because the ESS number-counts extend nearly 3 magni- 
tudes fainter than the spectroscopic catalogue, comparison 
of the observed galaxy counts with those predicted by the 
LFs provides a test of how well the measured LFs and 
evolution rates can be extrapolated from z ~ 0.5 to z ~ 1. 

Figures El ^2 El show the observed number counts 
in the R r-, V and B bands r esp., binned in intervals of 
0.25 ma s ■lArnonisei al .l ll997l) . In each band, the magni- 
tude of the faintest plotted point corresponds to the com- 
pleteness limit: R c ~ 23.5, V ~ 23.5, and B ~ 24, respec- 
tively. At bright magnitudes, the plotted number count 
distributions start in the first 0.25 mag bin where the count 
is larger or equal to 50 galaxies, as the counts are highly 
uncertain at low count level. For the same reason, the fits 
described in the following start in the first bin where the 
count is larger or equal to 100, that is in the bins centered 
at R c = 18.875, V = 19.625 and B = 20.325, respectively. 

For modeling the galaxy number-counts, one must de- 
fine a set of K-corrections. We cannot use the ESS poly- 
nomials K-corrections described in Sect. OH as these are 
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Fig. 9. K-corrections in the R c bands (top panel), V 
(middle panel) , and B (bottom panel) derived from 
IColeman et al.l l|l98fl labeled "CWW extrap" in graph; 
see text for details) which are used for predicting the ESO- 
Sculptor galaxy number counts in Figs. EI| ^2 an d El 
Elliptical, Sbc and Magellanic Irregular types are plot- 
ted as dashed lines (running from top to bottom), and 
are used to model the ESO-Sculptor K-corrections for the 
early-type, intermediate-type, and late-type galaxies re- 
spectively; K-corrections for "CWW extrap" type Scd are 
shown as dotted lines, for comparison. The polynomial 
K-corrections calculated for the 3 ESO-Sculptor spectral 
classes and valid only for z < 0.6 are shown as heavy 
solid lines, with early-type, intermediate-type and late- 
type galaxies plotted from top to bottom. 



unconstrained at z > 0.6, whereas the number-counts to 
magnitudes ~ 24 are produced by galaxies out to z ~ 1.5. 
In replacement, we u se the K-corrections o btained from 
the optical spectra of IColeman et all <)l980h . which have 
been extrapolated in the UV and the IR by using the 
theor etical SEDs of the GISSEL library ([Chariot et alJ 
Il996|) for the 4 types Elliptical, Sbc, Scd, and Magellani c 
Irregular (see ISawicki et all Il997t lArnouts et al.l ll999fl . 
These K-corrections are shown in Fig. ^\ for the B, V 
and i? c bands (labeled "CWW extrap"). We choose to use 
for the ESS early-type, intermediate- type, and late-type, 
the K-corrections for the Elliptical, Sbc, and Magellanic 
Irregular types respectively. This choice is motivated by 
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Fig. 10. Comparison of the ESO-Sculptor magnitude 
number-counts in the R c band per 0.25 mag interval (open 
and filled circles) with the modeled counts defined as the 
sum of the predicted number counts for the 3 considered 
spectral classes: early- type (red dotted line), intermediate- 
type (green dashed line), late- type with no-evolution (blue 
dot-dashed line), late- type with the best-fit evolution rate 
(blue solid line). The Gaussian+Schechter composite lumi- 
nosity functions listed in Tableware used for the 3 classes. 
The amplitude tf>o of the Gaussian component of the lu- 
minosity function used for each galaxy class is identical 
to that used in Fig. ^| The black heavy dot-dashed line 
show the predicted total counts with no-evolution, and the 
black heavy solid line the adjusted total counts with evo- 
lution in the late- type galaxies (see text for details) . Only 
the data points with more than 100 galaxies per bin (filled 
circles) are used for the adjustment of the total evolving 
counts. The plotted error bars in the observed counts are 
estimated as yjn. The inset shows the predicted redshift 
distributions for the 3 spectral classes, with the same line 
coding as for the number-counts. 

the comparison for z < 0.6 of the "CWW extrap" K- 
corrections with the ESS K-corrections for the 3 spec- 
tral classes, also shown in Fig. |5J In the common redshift 
interval, the 2 sets of K-corrections are in good agree- 
ment, with the largest deviations occurring in the B band 
(~ 0.8 mag for the early-type, ~ 0.3 mag for the intermedi- 
ate type, and ~ 0.1 mag for the late-type galaxies), as it is 



Fig. 11. Same as Fig. 1101 for the ESS number-counts per 
0.25 mag interval in the V band. The amplitude </>o of the 
Gaussian component of the luminosity function used for 
each galaxy class is identical to that used in Fig. [21 

the most sensitive band to the template shape in the UV 
at the redshifts considered here. The consistent evolution 
rates obtained in the following from the number-counts 
in the R c , V and B bands a posteriori indicates that the 
CWW-extrap K-corrections do not introduce any severe 
systematic effects in the modeled number- counts. 

We then model the expected number-counts in each 
band as the sum of the predicted number-counts over the 
3 galaxy spectral classes T$: 

N(m) = V 3 Ts=1 N Ts (m), (23) 

where m is the apparent magnitude in either of the R c , 
V or B filters. For each class, the expected number count 
Nt s (m) is obtained by integrating the corresponding com- 
posite LF $ Ts (M) (listed in TableUJ over all redshifts < 
z < -Zhigii contributing to absolute magnitude M(m, z, Tg) 
at a fixed value of m (see Eq. 

N Ts (m) = / $t s [M(m, z, T s )} — dz; (24) 
Jo dz 

is the comoving volume element per unit solid angle 
defined as 

dV = c_ d L (z) 2 1 f25) 

dz h (i + z) 2 ^n m {i + zf + n A 1 ' 
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Fig. 12. Same as Fig. ^] for the ESS number-counts per 
0.25 mag interval in the B band. The amplitude 4>o of the 
Gaussian component of the luminosity function used for 
each galaxy class is is identical to that used in Fig. IT51 



The values of Zhigh are 1.5 in the R c band, and 2.0 in the 
B and V bands. Note that N(m) and N Ts (m) in Eas.l2"3l 
and |^] resp. are defined as number counts per unit solid 
angle and per unit magnitude interval. The curves plot- 
ted in Figs. US to are then multiplied by 0.25 mag and 
7r 2 /180 2 (to convert to number counts per 0.25 mag interval 
per square degree). 

In Figs. E3 EH and H3 the thin dotted, dashed, 
and dot-dashed lines correspond to the predicted 
counts Nt s ijn) for the early-type, intermediate-type, 
and late-type galaxies without evolution, and using 
(n m ,OA)=(0.3,0.7). The amplitudes <fo> and <jf of the com- 
posite LFs are those which match the integrals of the ob- 
served and predicted redshift distributions in the interval 
0.01 < z < 0.81; these values are listed in Table [21 and 
are also indicated in Figs . 1 1 31 1 1 41 and Hoi fsee Sect.EJ. The 
summed number-counts over the 3 classes are plotted as 
heavy dot-dashed lines in Figs. EH and The to- 
tal non-evolving counts only match the observed number- 
counts at magnitudes brighter than R c ~ 21.5, V ~ 21.5 
and B ~ 22.0 resp., which corresponds or is close to the 
magnitude limit of the respective redshift samples. At the 
faint limit, the no-evolution counts under-predict the ob- 



served counts by ~ 25% in the R c band, and by ~ 60% in 
the V and B bands. 

Moreover, there is no value of non-evolving amplitude 
for the late-type LF which can both match the bright and 
faint ESS number-counts in all 3 bands. A simple scaling 
of the amplitude of the late-type LFs by a factor 1.37, 
1.85 and 1.90 in the P c , V and B bands resp., obtained 
by a weighted least-square minimization of the summed 
counts to the observed counts, still fails to match simul- 
taneously the observed number-counts at bright and faint 
magnitudes in all 3 bands: this scaling makes little change 
to the slope of the total counts, and essentially shifts them 
upward. 

Figures lTnillll andll2lshow that in each band, the num- 
ber counts at magnitudes fainter than the limit of the red- 
shift sample (~ 21 — 22 mag ) are dominated by the late-type 
galaxies. Introducing a scaling factor or some evolution 
in the amplitude of the early- type and/or intermediate- 
type LFs would therefore bring no improvement in match- 
ing the observed faint number-counts. Better adjustments 
may only be obtained by increasing the contribution from 
the late-type galaxies at faint magnitude. This provides 
further evidence that the late-type galaxies are evolving. 
In the following, we show that by introducing a linear (or 
power-law) evolution in the amplitude of the late- type LF, 
one obtains a very good adjustment of the observed num- 
ber counts in the 3 filters. Note that we deliberately do not 
consider any evolution in the early-type and intermediate- 
type ESS galaxy populations, for the following reasons: 

— there are several indications of a marked decrease 
in the number d e nsity of E galax ies with redshift 
(|Fried et alJ l200lt IWolf et all l2003|) which compen- 
sates for their passive luminosity evolution and ex- 
plains that no evolution is detected by the ESS in this 
population. 

— as far as the ESS intermediate spectral class is consid- 
ered, although a significant brightening due to passive 
evolution is expected for galaxies with present-day col- 
ors resembling those of Sb to Sbc morphological types 
(using intermediate v alues between th ose provided for 
Sa and Sc galaxies bv lPoggiantill99^ . little or no evo- 
lution in thei r luminosity density is detected out to 
z ~ 0.5 - 1.0 llLin et alJll993 IWolf et allEiol . 

We now determine the optimal evolution rate P0.15 
defined in Eq. |201 for the late-type galaxies by a 2-stage 
procedure. We first vary the value of P0.15 and per- 
form a weighted least-square fit of the summed predicted 
counts to the observed counts, with the amplitude of 
the early-type and intermediate-type LFs kept fixed. The 
weights are defined as the square-root of the observed to- 
tal counts, and might thus underestimates the true uncer- 
tainty in the observed counts, which should also account 
for galaxy clustering; we however verified that increasing 
the r.m.s. errors by as much as a factor of 2, a wide overes- 
timate of cosmic variance over the area of the ESS, would 
make negligeable change in the derived evolution rates. 
For each value of P0.15, the amplitude of the late- type LF 
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Table 3. Normalizing amplitudes 0o and <f>* of the Gaussian and Schechter components of the early-type, intermediate- 
type and non-evolving late-type luminosity functions of the ESO-Sculptor survey, in the Johnson-Cousins B, V and 
R c bands. 





(n m ,n A )=(o.3,o.7) 

Early-type Intermediate-type 

00 00 0* 


Late- type 

00 0* 


{n m ,n A )=( 

Early-type 

00 


1.0,0.0) 

Intermediate- type 

00 0* 


Late- type 

00 0* 


Rc < 21.5 
V < 21.0 
B < 22.0 


0.00310 0.00296 0.00387 
0.00323 0.00287 0.00375 
0.00344 0.00294 0.00385 


0.00339 0.03678 
0.00271 0.02939 
0.00299 0.03247 


0.00504 
0.00497 
0.00523 


0.00489 0.00639 
0.00444 0.00581 
0.00458 0.00600 


0.00520 0.05649 
0.00410 0.04450 
0.00457 0.04961 



Notes: 

- The listed values of 0o and 0* are in units of h 3 Mpc -3 mag -1 , and are obtained with Ho = lOO/i km s _1 Mpc" 1 . 

- 0o and 0* are calculated by normalizing the integral of the expected redshift distribution with the listed apparent magnitude 
limit to the observed number of galaxies in the interval 0.01 < z < 0.81 (see Sect.0. 

- The value of 0* for the intermediate-type and late-type galaxies is related to the corresponding value of 0o by the ratio 
0o/O.41nlO 0*, listed in Table0 



Table 4. Evolution rates and zero-points for the amplitudes 0o an d 4>* of the Gaussian and Schechter components of 
the late-type luminosity function in the ESO-Sculptor Survey, as derived from the number-counts in the Johnson P, 
V and Cousins R c bands. 







)=(0.3,0.7) 












,n A )=(i.o,o.o) 












-Po.15 


0o(O.15) 


0*(O.15) 


7 


00 (0) 


0*(O) 


Po.is 0o (0.15) 


0*(O.15) 


7 


00 (0) 


0*(O) 


Rc 


2.2 


0.00252 


0.02732 


1.8 


0.00211 


0.02292 


4.6 


0.00320 


0.03470 


2.5 


0.00287 


0.03119 


V 


4.3 


0.00216 


0.02347 


2.5 


0.00186 


0.02019 


5.7 


0.00322 


0.03493 


2.8 


0.00283 


0.03072 


B 


3.3 


0.00267 


0.02899 


2.2 


0.00226 


0.02449 


5.0 


0.00355 


0.03857 


2.6 


0.00313 


0.03399 



Note: 

- The listed values of 0o and 0* are in units of h 3 Mpc -3 mag -1 , and are obtained with Ho = lOO/i km s _1 Mpc 

- The linear, power-law parameterization of the evolution in 0o and 0* is defined in Eqs. 1201 1221 respectively. 

- Each value of 0* is related to the corresponding value of 0o by the ratio 0o/O.4 In 10 0* =0.1 (see Tabic 0. 



is defined by matching the observed and predicted late- 
type redshift distribution in the interval 0.01 < z < 0.81 
(see next Sect.). The reference amplitude of the late-type 
LF is therefore a function of Pq.15- As this value may 
however not provide the optimal match between the pre- 
dicted and observed late- type counts, we allow in each 
least-square fit for a scaling factor to the amplitude of the 
late-type LF. The full procedure yields a first estimate 
of the evolution parameter Po.15 (defined as the value for 
which the reduced x 2 is smallest): 2.2, 4.3 and 3.3 in the 
P c , V, and B bands resp., with scaling factors 0.931, 1.005 
and 1.1287 for the late-type LF amplitude. The predicted 
number-counts for the late-type galaxies with the above 
quoted evolution rates and scaling factors are plotted as a 
light solid line in Figs. 1101 II ll and lT2l and the correspond- 
ing summed counts over the 3 galaxy types as a heavy 
solid line. 

We have also looked for other minima of the reduced 
X 2 by allowing for a scaling factor in the amplitude of the 
early-type and intermediate-type LFs, and searching for a 
common scaling factor for the 3 galaxy types. For values 
of Po.15 around the first minima listed above, we iterate 



over the values of the scaling factors for the 3 classes: 
the output scaling factors for the late-type galaxies are 
applied to both the early-type and intermediate-type LF 
amplitudes, and the new scaling factor for the late-type 
galaxies which minimizes the x 2 is calculated. After 5 to 
10 iterations, this converges to a common scaling factor for 
the 3 galaxy types. When considering the final reduced % 2 
obtained by these iterations, a slightly smaller evolution 
rate Po.15 = 2.0 is obtained for the P c counts, with a com- 
mon scaling factor of 0.972; the same minimum Po.15 = 4.3 
is confirmed in the V band, with a common scaling factor 
of 1.003 for the 3 galaxy types; and a slightly higher evo- 
lution rate Po.15 = 3.8 is obtained in the B band, with a 
common scaling factor 1.055. Note that in each filter, the 
common scaling factor for the 3 galaxy types is closer to 
unity than the factor used when scaling only the late-type 
LF. The predicted number-counts with a common scaling 
factor are indistinguishable from those in which only the 
late- type LFs are scaled (which are shown in Figs. ll0ll"2*|l . 

The excellent adjustment of the number counts us- 
ing the linear evolution model of the amplitude of the 
late-type LF while keeping nearly constant the density 
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of early-type and intermediate-type galaxies provides ev- 
idence that the late- type galaxies evolve out to z ~ 1.0. 
The inserts of Figs. lTUlllll andll2lshow that in the 3 bands, 
the expected redshift distribution of the galaxies detected 
in the ESS number counts have a peak near z = 0.5 and 
extend to z > 1. Note also that the number counts provide 
better agreement among the evolution rates in the 3 bands 
(Po.is = 2.0 - 2.2 in P c , P .is = 4.3 in V, P„.i5 = 3.3 - 3.8 
in B) that those measured from the redshift survey only 
(Po.is = 3.5 in i? c , Po.is = 7.5 in V, P .15 = 8.6 in B; see 
Sect. El . 

We estimate the uncertainty in the value of P0.15 mea- 
sured from the number-counts by applying the following 
tests: (i) in the late-type LF, we change alternatively the 
slope a of the Schechter component from —0.3 to 0.39, the 
value actually measured from the P c < 21.5 sample (see 
Paper I), and the peak magnitude Mo of the Gaussian 
component by ±0.3 mag , as these are the 2 parameters 
which have the largest impact on P0.15; (ii) we vary the 
amplitude of the early-type or intermediate-type LFs by 
±10%, which provides a conservative estimate of the un- 
certainties in the values of (jf and 4>q for these samples 
(see Eqs. I16I17|) . Each of these tests yields a change in 
P0.15 by ^ 0.5. We thus adopt as a conservative uncer- 
tainty cr(Po.i5) ~ 1.0. The above values of P0.15 obtained 
from the number-counts then differ by 1 — 2a from filter 
to filter. 

We emphasize that despite the incompleteness in blue 
galaxies in the ESS V and B spectroscopic samples, when 
used together with the ESS V and B magnitude number- 
counts, they do provide useful constraints on the evolution 
rate for the late-type galaxies, and yield consistent results 
with those derived from the R c number-counts. Note that 
in the estimation of P0.15 from the number-counts, the 

V and B spectroscopic samples are used only to derive 
the amplitude of the LFs by normalizing to the observed 
redshift distributions. The consistent evolution rates ob- 
tained in the B, V, and R c bands reinforces the detected 
evolution as a real effect. The tendency of an increased 
evolution rate measured in the V and B bands compared 
to the P c , may be due to the higher sensitivity of the 

V and B bands to the late-type galaxies, already men- 
tioned in Sect. 15.31 at the peak redshift probed by the 
ESS number-counts (z ~ 0.5), taken together the B and V 
bands probe the rest-wavelength interval ~ 2700 — 4000A, 
lying just blue-ward of the Call H & K break; the signif- 
icant star formation activity present in late-type galaxies 
does cause an increased flux at these wavelengths. 

Note that applying the above analysis using pure 
Schechter LFs with a common slope a = —1.48 for the P c , 
V, and B LFs (Paper I) yields values of Po.is smaller by 
0.5 than those derived from the composite LFs. However, 
the use of pure Schechter LFs yields a marked degeneracy 
between the slope a of the late-type LF and the evolution 
rate Po.is- For example, changing the slope a of the late- 
type LF from —1.64, measured from the P c < 20.5 sample, 
to —1.48, measured from the R c < 21.5 sample (Paper I) 
yields an increase in P0.15 by nearly one unit. This is to be 



contrasted with the change of only 0.3 in P0.15 obtained 
when changing the slope of the Schechter component of 
the late-type composite LF from —0.3 to 0.39 (note the 
large variation). The degeneracy in the faint-end slope a 
of a Schechter LF could be partially reduced using the red- 
shift distribution, but there remains large uncertainties, as 
the ESS spectroscopic sample is far from a fair sample of 
Universe, and the redshift distribution fails in averaging 
out the large-scale structure (see Sect. 01 . 

In Table 01 we list the values of P0.15 obtained in the 
3 filters for (fi m ,f2A)=(0.3,0.7); there, the secondary iter- 
ation stage aimed at obtaining a common scaling factor 
for the 3 spectral classes is not used, as it makes a small 
difference in the evolution parameter. We also measure 
P0.15 f° r (^m,^A)=(l-0,0.0); in that case, the LF char- 
acteristic magnitudes Mq and M* listed in Table are 
shifted by AM = 0.3 mag , the variation in absolute magni- 
tude corresponding to the change in luminosity distance 
at z = 0.3 (the approximate peak redshift of the ESS; see 
Sect. UJ|. The corresponding values of <fio an d 4>* for the 
early-type, intermediate-type and non-evolving late-type 
LFs, obtained by normalizing the integral of the expected 
redshift distribution with (O m ,f2A)=(l-0,0.0) to the ob- 
served number of galaxies in the interval 0.01 < z < 0.81 
are also listed in Table 01 

Finally, we apply the power-law evolution model of 
Eg. 1221 to the number-counts and derive the best-fit value 
of 7 for both sets of cosmological parameters; the resulting 
values of 7 are listed in Tabled As for the P0.15 parameter, 
the uncertainty in 7 is estimated to be of order of 1.0. In 
the 3 bands, the minimum \ 2 is systematically smaller for 
the linear evolution model than for the power-law model, 
but the difference is small. Note that the larger evolu- 
tion parameters obtained with cosmological parameters 
(fi m ,QA)=(l-0,0.0) are due to the corresponding smaller 
volume element at increasing z. 

7. The ESS redshift distributions per 
spectral-type 

In Figs. 1131 1141 and 1 151 we compare the observed redshift 
distributions for the 3 spectral classes in the P c < 21.5, 
V < 21.0 and B < 22.0 samples resp. with the expected 
distributions calculated using the composite LFs listed in 
Table ^ F° r the observed distributions, we plot the 2 
histograms obtained with a redshift bin Az = 0.08 and 
offset by 0.04 in redshift, in order to illustrate visually 
the uncertainties in the observed distribution. A large 
bin size in redshift is used in order to smooth out the 
variations due to large-scale clustering; this scale would 
correspond to ~ 120 hr 1 Mpc at small redshift, larger 
than the typical size of the voids in the redshift surveys to 
z < 0.1 ide Lamwent et alJlTflSri ISnectman et a.ljll996t 
ISmall et alll997aHColless et al.l200lUZehavi et all2002h . 
and comparable to the scale of the largest inhomogeneities 
detected so far in r edshift surveys IjBroadhurst et al EM 
Ifteller et a.lj Ffl97L The marked deviations between the 
2 histograms are due to large-scale structure on even 
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Fig. 13. The observed redshift distribution for the early- 
type (top graph), intermediate- type (middle graph), and 
late-type galaxies (bottom graph) in the ESS R c < 21.5 
sample. In order to allow for the noise in the observed his- 
tograms, these are plotted for 2 binnings of Az = 0.08, off- 
set by 0.04 in redshift (first bin starting at z = 0, z = 0.04 
for the dotted, resp. solid histogram) . The expected distri- 
bution for a uniform galaxy distribution with the ESS in- 
completeness and the composite luminosity function listed 
in Table ^ is over-plotted as a solid line for the 3 spec- 
tral classes. For the late-type galaxies, we also plot as a 
dashed line the expected distribution with an evolution 
factor Po.15 = 2.2, as derived from the R c number counts 
(see Sect. El . All expected curves are normalized so that 
they match the observed number of galaxies in the inter- 
val 0.01 < z < 0.81; the resulting amplitudes <f>o and (jf 
for the Gaussian and Schechter components are indicated 
inside each panel. 



larger scales: a deficit of observed galaxies in the in- 
terval 0.33 z ^ 0.39, and an excess in the interval 
0.39 Si z 0.46; it is however unclear whether these struc- 
tures extend beyond the limited angular scale of the ESS. 

As in Sect. 15. ll the expected curves in Figs.ElEl and 
1X31 are based on the integral of the selection function over 
the bin-size Az = 0.08, and we use the K-corrections cal- 
culated for the average spectral-type < T§ > among each 
spectral class (see Eq. ^| . For the 3 galaxy types and in 
the 3 filters, the amplitudes 4>o an d 4>* defining each non- 
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Fig. 14. Same as Fig. US for the ESS V < 21.0 sample, 
with .Po.15 = 4.3 as derived from the V number counts 
(sec Sect. El. 



evolving expected curve are defined by normalizing the 
integral of the expected distribution to the observed num- 
ber of galaxies in the interval 0.01 < z < 0.81 (the ratio 
0o/O.41nlO (j)* takes the values listed in Table E}. The re- 
sulting values of 4>o and </>* are indicated inside each graph 
of Figs. El to El an d are also listed in Table El in column 
labeled (Q m ,Q\)— (0.3,0.7). These values of </>o differ from 
the estimates $3(0.51) listed in Table0 because the latter 
result from an integral over the narrower redshift interval 
0.1 < z < 0.51 (see Sect.ETTtf. 

For the early-type and intermediate-type galaxies, the 
expected distributions in Figs. El El an d El provide a 
good match to the observed histograms. For the late-type 
galaxies, the expected no-evolution redshift curves (with 
-Po.15 = 0.0) show a systematic shift towards low redshifts 
when compared to the observed distributions. Moreover, 
the expected curves lie systematically near the lower val- 
ues of the observed histograms for z <; 0.3. These effects 
are present in the 3 filters. 

For the late-type galaxies, we also plot in Figs. El 
El an( l El the expected redshift distributions with the 
values of the evolution factor Po.15 listed in Table 01 for 
(fi m ,Q A )=(0.3,0.7). The values of </> (0.15) and 0*(O.15) 
which normalize the integral of each evolving distribution 
to the corresponding observed number of galaxies in the 
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Fig. 15. Same as Fig. [H for the ESS B < 22.0 sample, 
with Po.is — 3.3 as derived from the B number counts 
(see Sect.EJ). 

interval 0.01 < z < 0.81 are indicated in Figs. 1131 to 1151 
Note that these values differ from those listed in Table 
as the latter are derived by normalization to the total 
number-counts. The difference is however small, 10%, 
thus bringing a posteriori evidence of consistency between 
the redshift and magnitude distributions. Note also that in 
the evolving curves, we have extrapolated to0<2:<0.15 
the linear evolution of $ for the late-type galaxies param- 
eterized in Eq. 1201 although it was measured from the re- 
stricted redshift interval 0.15 < z < 0.5. Our motivations 
for this choice are: 

— if we assume that $ remains constant and at the value 
of the linear fit at z — 0.15 for redshifts < 0.15, the 
expected redshift distributions show an marked excess 
of galaxies at these redshifts, which does not match 
the observed distribution (this effect is observed in all 
filters). 

— the location of the ESS was visually selected by ex- 
amining copies of the ESO/SERC R and J atlas sky 
survey with the criteria to avoid nearby galaxies and 
nearby groups and clusters of galaxies. This implies 
that the ESS survey has a systematically low density 
of galaxies at z Si 0.1. 

In contrast to the no-evolution curves, the evolving 
late-type distributions provide a good match to the ob- 
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Fig. 16. The observed redshift distribution for the bright 
and faint galaxies in the ESO-Sculptor late-type class 
(M(R C ) < -19.0 + 51og/i in top panel, M{R C ) > -19.0 + 
51og/i in bottom panel). In order to allow for the noise 
in the observed histograms, they are plotted for 2 bin- 
nings of Az = 0.08, offset by 0.04 in redshift (first bin 
starting at z = 0, z = 0.04 for the dotted, resp. solid his- 
togram). The expected distribution for a uniform galaxy 
distribution with the ESS incompleteness and the com- 
posite late-type LF listed in Table and truncated at 
M(i? c ) = —19.0 + 51og/i is plotted as a solid line in both 
panels. The expected distribution with an evolution fac- 
tor Pq.i5 = 3.0 is plotted as a dashed line in both panels; 
the expected curve with P0.15 = 6.0 is also shown as a 
dash-dotted line in the top panel. The expected curves 
are normalized so that they match the observed number 
of galaxies in the interval 0.01 < z < 0.81. 



served histograms in Figs. 1131 1141 and H"5l the low redshift 
peak is shifted to higher redshift (z *Z 0.2), and the high 
redshift tail has a higher amplitude which better matches 
the observed data. The best agreement of the observed 
and evolving expected curve is obtained in the R c band 
(Fig.[T3J|. In the V and B bands (Figs. [H and O, the ex- 
pected curves with P0.15 = 4.3 and Po.is = 3.3 resp. may 
be still be too low at z 0.3. Using Po.15 = 7.0 and 
Po.is = 8.0 in the V and B bands (as obtained from the 
direct fits of <f>i in Sect. 15. 3|) yields a better match of 
the observed redshift distributions at z ^ 0.3, but sys- 
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tematically underestimate the distributions at z ^ 0.15. 
This may be an indication that extrapolation of the lin- 
ear model of Eq. [20|to z < 0.15 is not satisfying for large 
values of Pq.is- The varying incompleteness with galaxy 
type which affects the V and B samples is also likely to 
complicate the adjustment of the redshift distributions. 

Note that the incompleteness of the full sample is 
used to correct for the incompleteness in apparent mag- 
nitude of each spectral-type sample in Figs. 1131 1141 and 
1151 Ideally, one should use the incompleteness calculated 
for each spectral class. However, galaxies with no red- 
shift measurement have no spectral-type determination. 
We have examined the dependence of the incompleteness 
as a function of galaxy-type using the colors of the galax- 
ies without redshift, as these are correlated with galaxy 
type. For the R c < 20.5 sample, the incompleteness is 
uniform with galaxy colors, justifying the use of the aver- 
age incompleteness for the full sample. For the V < 21.0 
and B < 22.0 samples, the incompleteness is significantly 
stronger for the bluer galaxies; we cannot however evalu- 
ate the incompleteness per spectral-type from the colors 
as the relation between color and spectral-type suffers a 
large dispersion. The relative larger incompleteness in blue 
galaxies at faint magnitudes of the V and B spectroscopic 
samples converts into a relative larger incompleteness in 
faint late-type galaxies compared to early-type galaxies. 
This might explain why the expected curves for the late- 
type galaxies in Figs. 1141 and 1151 appear to systematically 
under-estimate the observed distribution. 

As already mentioned in Sect. the ESS late-type 
class contains predominantly Sc+Sd and dl galaxies. An 
interesting issue is how the giant (Sc+Sd) and dwarf (dl) 
galaxies contribute to this evolution. In the composite fit 
to the late- type LF at R c < 21.5 shown in the lower left 
panel of Fig. ^ the giant and dwarf components cross- 
over near M{IQ ~ -18.5+5 log h for (O m ,O A )=(1.0,0.0); 
for (fi m , fiA)=(0.3,0.7), this magnitude converts approx- 
imately to M(R C ) ~ -18.8 + 5 log ft (see Sect. gj. We 
therefore calculate the observed and expected redshift 
distributions for the 2 following sub-samples: the bright 
sample, with M(R C ) < —19.0 + 5 log ft., containing 227 
galaxies, which are predominantly Sc+Sd galaxies; the 
faint sample, with M(R C ) > —19.0 + 5 log ft, containing 
82 galaxies, which are predominantly dl galaxies. Lower 
panel of Fig. ^] shows that an evolution factor with 
Po.i5 = 3.0 has a small impact on the distribution of 
galaxies with M(R C ) > —19.0+ 5 log ft, and both ex- 
pected curves are compatible with the observed distribu- 
tion. In contrast, the expected no evolution curve provides 
a poor match to the observed distribution of galaxies with 
M(R C ) < —19.0 + 5 log ft, whereas the expected curve with 
-Po.15 = 3.0 provides a better agreement. The expected dis- 
tribution for a higher evolution rate, Po.15 = 6-0, is also 
shown in Fig. llSl for the bright galaxies: it provides an even 
better adjustment (note that if the dl galaxies where not 
evolving, a higher evolution rate than measured for the 
Sc+Sd+dl altogether would be expected for the Sc+Sd 
galaxies alone). This suggests that the Sc+Sd galaxies are 



likely to contribute significantly to the detected late-type 
evolution. We however cannot exclude evolution of the dl 
galaxies. 

8. Comparison of the ESS evolution with other 
surveys 

Most other redshift surveys to z ~ 0.5 — 1.0 detect evo- 
lution in the luminosity function. The evolution affects 
either the amplitude <£> of the LF (number-density evolu- 
tion), or its shape via the characteristic magnitude (lumi- 
nosity evolution) and/or its faint-end behavior. Evolution 
in the luminosity density is also often used for measur- 
ing the evolution rate, as it has the advantage to account 
for the 3 types of evolution. Note that for a non-evolving 
shape of the luminosity function ip(M)dM , any evolution 
in the amplitude $ (see Eq.Ql yields an identical evolution 
rate of the luminosity density. 

In the following, we only consider the evidence for 
separate evolution of the E/S0 and Spiral galaxies, as 
evolution in the full galaxy population does not allow 
one to isolate the evolving population. For instance, the 
total luminosity density pj in the CFRS sh ows an in- 
crease with redshift which iLillv et all l)l99fifl model as 
p L {z) oc (1 + z) 2 7±0 - 5 at 4400 A for an (fi m ,n A )=(1.0,0.0) 
cosmology. This is close to the value 7 = 2.6 obtained from 
the ESS B band number-counts in Tabled suggesting that 
the evolution rate for the Spiral galaxies in the CFRS may 
be higher than in the ESS. A direct measure of the evolu- 
tion rate of the blue galaxies in the CFRS (bluer than a 
non-evolving Sbc galaxy) would however be required for a 
quantitative comparison with the ESS. 

Because there is weak evidence for evolution in 
the faint-end slo pe of the LF in the existing surveys 
iHevl et al.lll997h . we restrict the following discussion to 
the evidence for (i) number-density evolution, and (ii) lu- 
minosity evolution. 

8.1. Number-density evolution 

Ellis et al.l l|l996l) detect a marked density evolution in 
the Autofib star-forming galaxies by a factor 2 between 
z ~ 0.15 and z ~ 0.4, which would correspond to an over- 
all fading of the population by 0.5 mag in the bj band. This 
is comparable to the ESS evolution rate derived above: 
for P0.15 = 3.0, the density increases by a factor 1.75 
between z ~ 0.15 and z ~ 0.4. Further ana lysis of the 
Auto fib survey based on galaxy spectral types ( Hevl et all 
1997) leads to detection of a strong density evolution in 
the late-type spiral galaxies (Sbc and Scd) which they 
model as (f>*(z) oc (1 + z) 2 - 9 for the Sbc galaxies and as 



(j>*(z) tx (1- 



\3.5 



for the Scd galaxies in a cosmology with 



(fi m ,f2A)=(l-0,0.0). These values are in acceptable agree- 
ment with the 7 = 2.6 value derived from the ESS in the 
B band. 

In the CNOC2, the dominant evolution detected by 
iLin et al . (1999) in the redshift range 0.1 ^ z 0.55 is a 
strong density evolution of the late- type galaxies, defined 
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as galaxies with UBVR C I C colors similar to those com- 
puted from the Scd and Im templates of IColeman et alJ 
jl980() : the authors model the evolution as <fi*(z) oc e Pz , 
with P = 2.8 ± 1.0 and P = 2.3 ± 1.0 in the B AB band for 
fl m = 1.0, ft m = 0.2 resp.; and P = 2.9±1.0, P = 2.6±1.0 
in the R c band for Sl m = 1.0, TO = 0.2 respectively. At 
z Si 0.1, e Pz and (1 + z) 7 can be expanded into 1 + Pz and 
1 + 7z resp., which would yield a reasonable agreement at 
low redshift between the CNOC2 and the ESS in both the 
B and R c bands. However, at z ~ 1, e 2 5z is a factor of 2 
larger than (1 + z) 2 - 5 , implying a stronger evolution rate 
in the CNOC2 than in the ESS. 

Obtained with a similar observational technique as the 
CNOC2, the field sample of the CNOC1 llLin et al. II 19971) 
shows an increase of the luminosity density of galaxies by 
a factor 3 between z ~ 0.25 and z ~ 0.55 for galaxies with 
rest-frame colors bluer than a non-evolving Sbc galaxy. In 
the NORRI S survey of the Corona Borealis Supercluster, 
ISmall et"aT1 l|l997bh detect a similar evolution rate: the 
amplitude 0* of the Schechter LF for galaxies with strong 
[Oil] emission line increases by nearly a factor 3 from z Ss 
0.2 to 0.2 < z < 0.5. These various values of the number- 
density evolution rate are stronger than the increase by a 
factor 1.8 in the ESS density of late- type galaxies which 
we derive for Po.is = 3.0 between z ~ 0.2 et z ~ 0.5. 

In the CADIS survey, in which redshifts are derived 
from a combinat ion of wide and medium-band filters, 
iFried et al.l l|200lj) detect an increase of the Johnson B 
luminosity density of Sa-Sc galaxies with redshift, which 
is partly due to an increase in the amplitude (jf of the 
fitted Schechter luminosity function, and can be mod- 
eled as p(z) oc (1 + 6.9z) for (Q ro ,f2 A )=(1.0,0.0), and 
as p(z) oc (1 + 0.78z) for (Q m ,fJ A )=(0.3,0.7). The evolv- 
ing term in the luminosity density can be converted into 
p(z) oc [1 + 3.4(2-0.15)], resp. p{z) oc [1 + 0.70(z- 0.15)]. 
Whereas the CADIS evolution rate is similar to that in the 
ESS for (f2 m ,f2 A )=(l. 0,0.0), it is much smaller than in the 
ESS for (f2 ro ,fi A )=(0.3,0.7) (see Table [P; note however 
that this comparison may be complicated by the fact that 
the considered evolving CADIS population contains early- 
type Spiral galaxies, to the contrary of the ESS late-type 
spectral class. 

In contrast, the recent COMBO- 17 survey l|Wolf et alJ 
2003) which is also based on a combination of wide and 
medium-band filters, shows no significant evolution in the 
number-density and luminosity density of either Sa-Sc 
and Sbc-Starburst galaxies from z ~ 0.3 to z ~ !•! i n 
the Johnson B and SDSS r bands (|Fukugita et alJll996l) . 
Moreover, whereas the various mentioned surveys (CFRS, 
Autofib, NORRIS, CNOC2, CNOC1, CADIS) show no 
or a weak change in the luminosity density of early-type 
(E-S0) or red galaxies over the considered redshift range, 
the COMBO-17 survey detects a marked increase with 
redshift by a factor of 4 in the contribution from the 
E-Sa galaxies to th e r and B lumi nosity densities for 
(n ro ,Q A )=(0.3,0.7) llWolf et allEool . The different re- 
sults between the COMBO-17 and the other redshift sur- 
veys may be due to the complex selection effects inher- 



ent to surveys based on multi-medium-band photometry 
such as the COMBO-17, and which are most critical for 
emission-line galaxies. These effects however do not seem 
to affect the CADIS survey. 

At last, the evolution detected in the far infrared from 
IRAS galaxies JSaunders et alJll990t iBertin et all Il997t 
iTakeuchi et all2003j) can be characterized as (j>*(z) oc (1 + 
z) 7 with 7 ~ 3 — 3.4 for pure density evolution (both 
cosmologies considered in this article are used, depending 
on the authors). The evolution of the IRAS galaxies is 
consistent with the ESS late-type evolution, in agreement 
with the fact that IRAS galaxies may represent a sub- 
population of the optical spiral galaxies. 

8.2. Luminosity evolution 

The apparent density evolution detected in the ESS could 
also be produced by a luminosity evolution of the late-type 
spiral galaxies: if these galaxies were brighter at higher 
redshift, they would enter the survey in larger numbers 
at a given apparent magnitude. Using the values of the 
power-law index 7 listed in Table |H an d the slopes o f the 
ESS magnitude number-counts l|Arnouts et aljfl997|) . we 
can estimate a corresponding magnitude brightening. For 
cosmological parameters (fi m ,fi A )=(l. 0,0.0), we measure 
from Tableman increase in the number-density of galaxies 
by a factor 2.8 between z ~ and z ~ 0.5, and by a factor 
5.7 between z ~ 0.5 and z ~ 1 in the R c band; in the B 
bands, the density increases by a factor 2.9 at z ~ 0.5, 
and 6.1 at z ot 1. For (fi m ,fi A )=(0.3,0.7), the density in- 
creases by 2.1 and 3.5 at z ~ 0.5 and z ~ 1 resp. in R c , 
and by 2.4 and 4.6 at z ~ 0.5 and z — 1 resp. in B. Using 
the slopes (3 = 0.38 in the R c band and (3 — 0.46 in the 
B band (n(m) oc 10^ m ) for the ESS magnitude number- 
counts l(Arnouts et alJll997T) . these values of the density 
increase are equivalent to ~ 1.2 mag and ~ 2.0 mag bright- 
ening of the late-type ESS galaxies at z ~ 0.5 and z ~ 1 
resp. in R c , and a ~ l.i ma s and a ~ 1.7 mag brightening 
resp. in B for (n m ,ft A )=(1.0,0.G); and to - 0.8 mag and 
~ 1.4 mag brightening resp. in both the R c and B bands 
for (ft m ,fi A )=(0.3,0.7) (the steeper evolution rate in the 
B band is compensated by a steeper slope of the number- 
counts also in the B band). 

These brightening estimates for the ESS late-type 
galaxies with (f2 m ,f2A)=(0.3,0.7) are comparable to those 
caused by the passive evolution of an Sc galaxy (due to 
the evolution of the stel l ar po pulation). From the model 
predictions of IPoggiantil (|l997L using (fi m ,fi A )=(0. 45,0.0) 
and Hq = 50 km s _1 Mpc -1 , which imply an age of 
the Universe of 15 Gyr), an Sc galaxy brightens by ~ 
0.6 mag and - 1.2 mag in its rest-frame R band at z ~ 0.5 
and z ~ 1.0 resp., and by ~ 0.9 mag , ~ 1.5 mag resp. in 
its rest-frame B band. A comparable or stronger bright- 
ening is expected for the Elliptical galaxies at these red- 
shifts: ~ 0.6 ma g and ~ 1.4 ma g in R, ~ 0.7 mag and 
- 3.0 mag in B l|Poggiantilll997h . However, there are indi- 
cations of a marked decrease in the number density of E 
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galaxies with redshift l|Fried et al.l200lHWolf et alJl2003|) . 
which compensates for their luminosity evolution, which 
explains why no evolution in this population is detected 
by the ESS. 

iLillv et all l)l995l) detect a l mag brightening of the 
CFRS blue galaxies (defined as galaxies with rest-frame 
colors bluer than a non-evolving Sbc template from 
IColeman et al.lfl980h between the intervals 0.2 ^ z ^ 0.5 
and 0.5 ^ z ^ 0.75 in an (f2 m ,f2A)=(l. 0,0.0) cosmology; 
they however cannot discriminate whether this brighten- 
ing is due to luminosity or density evolution. With 7 = 2.6 
(see Table 0J, the ESS B density of late- type galaxies in- 
creases by a factor 1.6 between z ~ 0.25 and z ~ 0.625 
(the median value of the 2 quoted CFRS intervals), which 
corresponds to a brightening of ~ 0.44 mag , nearly a factor 
2 smaller than in the CFRS. ICohenI l)2002|) also detect for 
the emission-line dominated g alaxies at z ~ 1, when com- 
pared to the measurement of iLin et al.l l|l996|) at z ~ 0, 
a mild brightening by ~ 0.75 mag in the R band for an 
(r2 m ,f2 A ) = (0.3,0.0) cosmology, which is a factor 2 smaller 
than in the ESS estimated brightening in the R c band for 
(n m ,O A )=(0.3,0.7). 

9. Conclusions, further discussion and prospects 

Using the Gaussian+Schechter composite LFs measured 
for the ESO-Sculptor Survey, we obtain evidence for evo- 
lution in the late spectral-type population containing late- 
type Spiral (Sc+Sd) and dwarf Irregular galaxies. This 
evolution is detected as an increase of the galaxy density 
n(z) which can be modeled as n(z) oc 1 + Pq,15(z — 0.15) 
with P0.15 ~ 3 ± 1 or as n(z) oc (1 + z) 7 with 7 ~ 
2 ± 1 using the currently favored cosmological parameters 
(n m ,O A )=(0.3,0.7); for (Q ro ,ft A )=(1.0,0.0), P .ib ~ 4± 1 
and 7 ~ 2.5 ± 1. Both models yield a good match of the 
ESS BVR C redshift distributions to 21 - 22 mag and the 
number-counts to 23 — 23.5 mag , which probe the galaxy 
distribution to redshifts z ~ 0.5 and z ~ 1.0 respectively. 
Using both the redshift distributions and the number- 
counts allows us to lift part of the degeneracies affecting 
faint galaxy number counts: the redshift distributions al- 
low to isolate the evolving populations, whereas the faint 
number counts provide better constraints on the evolution 
rate. These results are based on the hypothesis that the 
shape of the LF for the ESS late- type class does not evolve 
with redshift out to z ~ 1. 

Examination of the other existing redshift surveys to 
z ~ 0.5 — 1.0 indicates that a wide range of number-density 
evolution rates have been obtained. The evolution rate 
of the Sc+Sd+dl galaxies detected in the ESS is among 
the range of measured values, with some surveys having 
weaker of higher evolution rates. The most similar survey 
to the ESS, the CNOC2, yields a twice larger increase 
in the number density of late-type Spiral and Irregular- 
galaxies at z ~ 1. 

A priori, density evolution indicates that mergers could 
play a significant role in the evolution o f late- t ype Spiral 
and Irregular galaxies. iLe Fevre et al.l l|2000h detect a 



~ 20% increase in the fraction of galaxy mergers from 
z ~ to z ~ 1, which can be modeled as oc (1 + z) 3 - 2 ; 
interestingly, examination of their Fig.l indicates that a 
significant fraction of the merger galaxies have a Spiral or 
Irregular structure. 

The ESS density increase for the Sc+Sd+dl galaxies 
could also be caused by a ~ i ma s brightening of these 
galaxy populations at z ~ 0.5 and a ~ 1.5 — 2.0 mag bright- 
ening at z 1 (depending on the filter and cosmological 
parameters). This luminosity evolution is compatible with 
the expected passive brighteni n g of Sc galaxie s at increas- 
ing redshifts l)Poggiantilll997f ). [Driver] ll200ll) also sh ows 
that the Hubble Deep Field fWilliams et alJ Il996h bi- 
variate brightness distributions for Elliptical, Spiral, and 
Irregular galaxies are all consistent with passive luminos- 
ity evolution in the 3 redshift bins 0.3 — 0.6, 0.6 — 0.8, 
0.8 — 1.0. The ESS brightening at z ~ 1 agre es with the 
value measured from the CFRS blue galaxies llLilly et al. 
1995), but is twice smaller than that measured bv ICohenI 
(2002) for emission-line dominated galaxies. 

In all analyses of the redshift and magnitude dis- 
tributions, the major difficulty is to distinguish be- 
tween luminosity and density evolution, as these pro- 
duce the same net effect on the redshift and magni- 
tude distributions. Interpretation of density and lumi- 
nosity evolution of a galaxy population is also com- 
plicated by possible va riations in the st ar formation 
rate with cosmic time: ILillv et alJ l|l99^ evaluate an 
increase in the star formation rate of galaxies with 
large disks by a factor of ~ 3 at z ~ 0.7, which 
shows as an increase of the l uminosity density at bluer 
wave le ngths. Using PEGASE llFioc fc Rocca-Volmerangel 



Il997h . |R occa-Volmerange fc Fiod <ll999D also show that 
the Sa-Sbc galaxies have a star formation rate which varies 
more rapidly in the interval ^ z ^ 1 than for the E/S0 
or Sc-Im galaxies. 

The ESS suggests that the Sc+Sd galaxies are an 
evolving population, but evolution in the dl galaxies can- 
not be excluded. Whether the Spiral galaxies, or the 
Irregular galaxies, or both populations contribute signif- 
icantly to the excess number-counts is still unclear from 
the various exist ing analyses: using photometric redshifts, 
iLiu et ail l)l998|) detected a significant amplitude increase 
and brightening of the U LF of Sbc and bluer galaxies in 
the redshift interval z ^ 0.5, together with an excess 
population of starburst galaxies at z ^ 0.3 which are ab- 
sent at z & 0.3 ; galaxy number cou nts in the near infrared 
l|MartiniE00ll iTotani et aL1l200lL which arc less sensi- 
tive to current star formation, also allow some moderate 
number evolution of the Spiral galaxies and/or Irregular 
galaxies (with 7 ~ 1, see above) . Evolutionary effec ts are 
also detected at higher redshifts: iDriver et al.l l l998h show 
that at 1 ^ z ^ 3, the redshift distributions of Sabc galax- 
ies to I > 25 and of Sd/Irr galaxies to / > 2 4 require som e 
number and luminosity evolution (but see lDriverll200l|) . 
This could be consistent with the recent deep optical and 
infrared observations which favor mi ld luminosity gvolu- 
tion of the overall galaxy population l(Pozzetti et al.l l2003: 
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iKashikawa et al.ll2003j) . Nevertheless, none of these sur- 
veys allow one to isolate one evolving population among 
the Spiral and Irregular/Peculiar galaxies. 

In contrast, some surveys favor Irregular/Peculiar 
galaxies as a major contributor to the excess count: an 
excess population of Irregular/Peculiar galaxies was di- 
rectly identified as the cause for a strong deviation from 
no-evolution in the J and K number-counts per morpho - 
logical type l|Glazebrook et all 1199.4 iHuang et af]ll998ft : 
a population of gently-evolving starbursting dw arves was 
also invoked to explain these excess objects JCamposI 
Il997ft . Using morphology of galaxie s obtai ned from Hubble 
Space Telescope images, llm et all i|l999|) provide further 
evidence for a marked increase in the relative abun- 
dance of Irregular and Peculiar galaxies which they in- 
terpret as starbursting sub-L* E/SO and Spiral galaxies. 
iTotani fe Yoshiil l|l998j) also detect such an excess popu- 
lation at z ^ 0.5. Using again Hubble Space Telescope 
imaging for galaxy morphology. iBrinchmann et al.l 1 1998) 
estimate a ~ 30% increase in the proportion of galaxies 
with an irregular morphology at 0.7 j$ z ; $ 0.9 (see also 
lAbraham etaDll996t IVolonteri et"alll2000|) . 

Other surveys identify Spiral galaxies as the evolv- 
ing population: morphological number counts based 
on Hubble Space Telescope images indicate that the 
Spiral counts rise more ste eply than the no-evolution 
model ({Abraham et al.lll996[) . By complementing space- 
based images with groun d-based spectroscopic redshifts, 
IBrinchmann et aTl l)l998[) d etect a ~ l m ^s brig htening of 
the Spiral galaxies by z ~ l. lSchade et alJ ||199(t() also mea- 
sure a ~ 1.6 mag brightening of the central surface bright- 
ness of galaxy disks at 0.5 ^ z ^ 1.1. 

We emphasize that measuring a reliable evolution rate 
requires a realistic parameterization of the intrinsic LFs 
of each galaxy population. We show here that the ESS 
Gaussian+Schechter composite LFs provide more robust 
constraints on the evolution rate than pure Schechter LFs, 
as a small change in the faint-end slope has a large inci- 
dence onto the number-counts. In contrast, variations in 
the faint-end slope of the dE and dl galaxies, the only 
populations for which the Schechter faint-end is poorly de- 
termined have a smaller impact on the adjustment on the 
number counts. ITotani et al.1 l)200l|) also show that sepa- 
rating the E and dE galaxies into a Gaussian+Schechter 
LF does yield a better agreement of the K number-counts 
at K 22.5 than a pure Schechter function for the joint 
class of E+dE. 

By using better measures of the intrinsic LFs for 
the various galaxy types, one should be able to ob- 
tain improved measurements of the galaxy evolution 
rates, and eventually to use the f aint numbe r coun ts to 
probe the cosmological parameters. iKoo et all lfl993'l and 
iGronwall & Kool (|l995h had already suggested that pas- 
sive luminosity evolution is sufficient to match the B num- 
ber counts, thus implying that the numbe r counts could 
be use d to probe the curvature of space. IPozzetti et alJ 
l|l996h and iMetcalfe etaD lj200lh then showed that pas- 
sive luminosity evolution with a low value of Q, m provide 



good fits to UV-optical- near- infrared galaxy counts. The 
recent adjustme nts of very deep optica l and near-infrared 
galaxy counts llTotani & Yoshii l l2000t IN agashi ma et alJ 
1200 it ITotani et alJ l200ll INagashima et all l2002h which 
probe the galaxy distribution to z ~ 3 further confirm that 
the Einstein-de Sitter cosmology with a)=(1.0,0.0) 
is excluded at a high confidence level (see also lTotani et alJ 
Il997t iHe etaLll2000h : whereas by using the currently fa- 
vored values of fi m = 0.3 and Oa = 0.7, these fits con- 
strain the evolution rate of galaxies in the hierarchical 
clustering picture, with only s ome mild number evolution 
of Sbc/Sdm galaxies allowed l|Totani et alJl20oij) . 

Under-going deep redshift surveys to z <; 1 raise 
new prospects along this line. The present ESS analy- 
sis shows the usefulness of using both the magnitude 
and redshift distributions for studying galaxy evolution. 
By obtaining the redshift distributions per galaxy type 
to z <; 1 over large volumes which average out the 
large-scale struc ture, spectroscopic red shift survey such 
as the VIRMOS llLe Fevre & Vettolanill2003h and DEEP2 
llDavis et al . 2003) projects should provide improved clues 
on the evolving galaxy populations at z ~ 1 and bet- 
ter constrain the nature of this evolution. If these surveys 
confirm that all galaxy types only experience mild lumi- 
nosity/density evolution, complementing the deep infrared 
galaxy counts with morphological information might allow 
one to lift the degeneracy in the faint infrared number- 
counts and to confirm whether a low matter-density uni- 
verse is favored. This will however require a detailed 
knowledge (i) of the luminosity-size relation for distant 
galaxies of various morphological type (used to model the 
selection effects cause d by the cosmological dimming in 
surface brightness, see lTotani fe Yoshiill200o|) . and (ii) of 
the interstellar and intergalactic extinction. 
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